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operator  formalism,"  developed  previously  by  this  author,  is  an  essential 
ingredient  in  the  derivation  of  the  final  equations  presented  herein,  even 
though  these  equations  could  be  obtained  without  the  use  of  that  formalism. 

• The  final  equations  of  motion  are  independent  of  the  material  properties 
of  the  body  under  consideration.  In  order  to  make  specific  use  of  these 
equations,  specific  constitutive  equations  must  be  postulated.  The  most 
typical  constitutive  equations  used  in  the  spacecraft  dynamics  literature 
is  to  consider  the  deformable  body  to  consist  of  a collection  of  hinged 
rigid  or  linearly  elastic  sub-bodies;  then  the  deformation  within  one  of 
these  sub-bodies  is  small,  but  the  deformation  between  the  sub-bodies 
may  be  arbitrarily  large. 
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This  report  is  part  of  a continuing  effort  in  the  Controls 
Department  of  The  Aerospace  Corporation,  of  developing  tech- 
niques for  analyzing  and  simulating  the  dynamics  and  control  of 
large  flexible  spacecraft.  The  following  are  previous  related 
publications. 


• Jerkovsky,  W.  , "The  Transformation  Operator  Approach 
to  Multi-Body  Spacecraft  Dynamics.  Volume  I,  The 
Transformation  Operator  Formulation;  Part  1,  Momentum 
Approach,  " Aerospace  Corp.  , El  Segundo,  Calif.  , 
TOR-0075(5624-01)-1,  6 December  1974. 


• Jerkovsky,  W.  , "The  Transformation  Operator  Approach 
i ; to  Multi-Body  Dynamics,  " Aerospace  Corp.  , El  Segundo, 

Calif.,  TR-0076(6901-03)-5,  10  May  1976. 

I Additional  related  reports  are  in  preparation, 

j ■ The  present  report  describes  work  performed  under 

I Engineering  Methods  of  the  Guidance  and  Control  Division, 

Engineering  Science  Operations.  I wish  to  thank  Ms.  Karen  Saito 
' for  her  skillful  typing  of  the  manuscript. 
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I.  INTRODUCTION 


During  the  last  twenty  years,  the  control  system  analyst  has 
been  faced  with  modeling  the  dynamics  of  multibody  spacecraft 
with  ever  increasing  fidelity.  This  requirement  for  increased 
dynamics  modeling  sophistication  is  due  in  part  to  increasing 
control  accuracy  requirements,  increasing  spacecraft  maneuver- 
ing requirements,  increasing  spacecraft  size  and  structural 
flexibility,  and  finally,  an  increasing  number  of  moving  parts 
on  board  the  spacecraft.  A considerable  number  of  papers^ 

has  been  devoted  to  this  subjected  prior  to  1970,  and  an  even 
17-44 

larger  number  since  1970.  The  dynamics  equations  can 

get  so  complicated  (because  of  the  large  number  of  terms)  that 
it  is  difficult  to  see  the  forest  for  the  trees.  A similiar  prolifer- 
ation of  terms  in  dynamics  equations  has  occurred  in  the  study  of 

45-  54 

the  dynamics  of  mechanisms  and  linkages 

The  purpose  of  the  present  paper  is  to  describe  the  general 
form  of  the  dynamics  equations  of  motion  for  an  arbitrary  de- 
formable body  undergoing  large  deformations.  In  order  to  keep 

the  equations  completely  general,  no  particular  material  proper- 

55-56 

ties  (i.e.  no  constitutive  equations  )will  be  introduced  . In 

57 

this  respect,  the  present  paper  is  similiar  to  that  of  McDonough 

except  that  we  do  not  introduce  any  equations  from  continuum 

mechanics.  Instead,  we  consider  the  deformable  body  to  consist 

25 

of  a large  number  of  particles  (say  10  molecules)  of  negligible 
58 

inertia  . The  deformable  body  consisting  of  N particles  then 

has  6 ''external''  and  n "internal"  degrees  of  freedom,  where 
n ^ 3N  - 6;  if  there  are  no  constraints  among  the  N particles 
then  n equals  3N  - 6,  but  if  there  are  s constraints  then 
n = 3N  - 6 - s.  The  separation  into  "external"  or  "rigid  body" 
degrees  of  freedom  and  "internal"  or  "deformation"  degrees  of 
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freedom  is  similiar  to  that  introduced  by  Teixeira  and  Kane 
If  the  equations  of  motion  are  linearized  in  the  deformations,  then 
our  equations  reduce  to  the  corresponding  equations  in  the  litera- 
ture. If  the  deformations  are  taken  to  be  zero,  then  we  recover 
the  equations  of  motion  for  a rigid  body.  A topological  tree  of 
rigid  bodies  can  be  handled  by  letting  the  deformation  coordinates 
represent  the  relative  gimbal  angles  (and  relative  hinge  displace- 
ments, if  any)  between  contiguous  rigid  bodies.  Almost  all  the 
spacecraft  dynamics  equations  in  the  literature  can  be  obtained  by 
introducing  appropriate  specializations  into  the  equations  present- 
ed herein.  However,  there  are  two  types  of  spacecraft  dynamics 
problems  not  considered  herein:  the  deployment  of  roll-out 
booms  or  solar  arrays,  and  mass  expulsion  as  in  rocket  engine 
firing. 

There  are  a number  of  controversies  and  misconceptions  in 
the  spacecraft  dynamics  literature  on  several  basic  mechanics 
issues.  One  of  these  controversies  deals  with  whether  a 
"Newton-Euler " (or  "vectorial")  dynamics  formulation  is  prefer- 
able to  a "Lagrange-Hamilton"  (or  "analytical  dynamics"  or 
"variational  dynamics")  formulation.  Often  it  is  claimed  that  a 
"Newton-Euler"  formulation  has  the  advantage  of  being  "more 
physical",  but  a "Lagrange-Hamilton"  formulation  is  preferable 
because  constraint  forces  are  automatically  eliminated.  We  will 
show  in  this  paper  that  our  approach  (which  might  be  called 
"algebraical  dynamics")  has  both  of  the  above  advantages.  Our 
equations  have  the  form  of  those  of  Newton  and  Euler,  but  the 
procedure  we  use  to  get  the  equations  is  essentially  due  to 

Lagrange  and  Hamilton.  Our  approach  is  similiar  to  that  of 

6)  3 ” 6 5 

Kane  except  we  make  extensive  use  of  transformation 

operators^^  ^^to  make  velocity  transformations,  which  in  turn 
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induce  appropriate  momentum  and  force  transformations.  The 
transformation  operator  formalism  is  based  on  Kron's  method  of 
subspaces  which  in  turn  is  based  on  the  procedure  used 

by  Lagrange  to  go  from  Newton's  vector  equations  for  particles 
to  Lagrange's  equations  in  terms  of  generalized  coordinates. 

Another  controversy  centers  about  whether  it  is  preferable 

to  use  velocity  or  momentum  as  state  variables.  Russell  has 

been  the  principal  proponent  of  the  momentum  formulation  of 

74-76 

spacecraft  dynamics  problem.  Recently,  Vance  and  Sitchin 
have  also  been  advocating  the  momentum  formulation.  The  con- 
troversy here  is  similiar  to  the  question  of  whether  Lagrange's 
or  Hamilton's  equations  are  preferable  for  spacecraft  dynamics 
equations.  We  will  see  that  our  approach  allows  us  to  get  the 
Lagrange  and  the  Hamilton  equations  in  a simple  algebraic  man- 
ner which  does  not  require  forming  Lagrangians  and  Hamilton- 
ians and  then  taking  partial  derivatives.  Thus,  our  approach 
allows  us  to  get  both  the  velocity  and  the  momentum  equations 
for  the  deformable  body. 

Nonlinear  dynamics  equations  for  a deformable  body  can  be 
obtained  by  expressing  the  particle  displacements  (relative  to  a 
fixed  or  floating  reference  frame)  in  terms  of  the  eigenfunctions 
or  eigenvectors  of  the  corresponding  (or  closely  related)  linear- 
ized equations  of  motion.  Different  choices  of  reference 
frames^®’ and  different  choice  of  eigenfunctions  or 
eigenvectors  can  be  made.  We  will  maintain  generality  by  not 
specifying  the  form  or  nature  of  the  particle  displacement  field; 
we  merely  assume  that  some  displacement  field  function  does 
exist,  and  we  assume  this  function  is  continuous  and  has  con- 
tinuous partial  derivatives  of  up  to  the  second  order  (so  that  we 
can  interchange  the  order  of  two  partial  differentiations). 


1 


The  dynamics  equations  for  a spacecraft  are  sometimes 
written  relative  to  a point  fixed  in  the  structure,  and  sometimes 
they  are  written  relative  to  the  total  center  of  mass.  We  will 
maintain  generality  by  writing  them  relative  to  an  arbitrary 
point  which  may  be  fixed  in  the  structure  or  not;  in  particular, 
the  point  may  be  the  total  center  of  mass. 

Our  dynamics  equations  for  a deformable  body  are  given  in 
Section  V,  page  65.  Our  "momentum  formulation"  equations  are 
G + X = K and  G = P • ^,  where  ct  is  the  system  velocity,  G is  the 
system  momentum  (which  is  defined  so  that  the  kinetic  energy  is 
given  by  T = y G*'  'a),  P is  the  system  mass  (which  yields 
T = a • M:  • cr),  K is  the  system  force  (which  yields  T = K 'a) 
and  X is  an  extra  term  which  must  be  added  to  the  time  derivative 
of  G so  that  the  result  is  K.  Our  "velocity  formulation"  equation 
is  11  ' CT  + Y = K where  Y is  an  extra  term  which  must  be  added  to 
the  system  mass  times  the  system  acceleration  so  that  the  result 
is  K.  Of  course,  our  dynamics  equations  are  not  adequately  de- 
scribed until  each  of  the  variables  in  the  equations  is  well  defined. 
The  variables  for  the  "external"  or  "rigid  body"  degrees  of  free- 
dom are  described  in  Section  II  in  terms  of  the  corresponding 
variables  for  the  individual  particles  which  constitute  the  deform- 
able body.  Section  III  describes  the  variables  for  the  "internal" 
or  "deformation"  degrees  of  freedom.  Most  of  Section  IV  is 
taken  up  in  determining  first  the  general  forms,  and  then  the 
explicit  expressions  for  X and  Y;  the  expressions  for  p and  K 
are  actually  quite  straightforward. 
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II.  CLASSICAL  VECTOR  EQUATIONS 


We  will  start  out  with  Newton's  law  for  a particle,  written 
relative  to  an  inertial  origin.  This  law  applies  to  each  of  the 
particles  of  the  system  of  particles  which  represents  our  de- 
formable body.  By  assuming  that  action  is  equal  and  opposite  to 
reaction,  the  internal  forces  drop  out  of  the  equations  which 
describe  the  external  degrees  of  freedom;  however,  chese  inter- 
nal forces  do  enter  the  equations  for  the  internal  degrees  of 
freedom  (except  that  purely  constraint  forces  drop  out). 

A.  EQUATIONS  RELATIVE  TO  AN  INERTIAL  ORIGIN 

In  this  subsection,  we  introduce  the  vector  notation  and 
80-81 

equations  which  we  use  to  describe  the  deformable  body 

relative  to  an  inertial  origin.  We  number  the  particles  from 
1 to  N and  we  refer  to  a typical  particle  as  the  i''^  particle  or 

particle  i,  where  1 s i s N.  Let  r.  be  the  position  vector  of  the 

th  ^ , 

i*^”  particle  relative  to  an  inertial  reference  origin,  r.  will,  in 

th  ^ ^ 

general,  vary  with  time,t.  The  velocity  of  the  i particle,  v^^, 
is  equal  to  the  inertial  time  derivative  of  and  is  denoted  by 


V.  = r. 

1 1 


This  is  a "material"  time  derivative^^  since  it  is  the  velo- 
city experienced  by  a particular  material  point  of  the  deforma- 
ble body. 

th 

Let  m^^  be  the  mass  of  the  i particle  and  let  M be  the 
total  mass  of  all  particles.  Then  a point  c,  the  total  center  of 
mass,  is  given  by  the  position  vector  where 
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j .-L  E 


M i 


where 


= E 


In  Equations  (2)  and  (3)"'it  is  understood  that  the  summation  is 
over  all  the  particles;  i.e.  for  i = 1 to  N.  The  velocity  of  the 


center  of  mass,  v^,  is  given  by 


V = r = — m.  V. 

^ = M i 


The  linear  momentum  of  particle  i,  and  of  the  whole 

-> 

system,  P,  are  given  by 


p.  = m.  V. 

11 


Combining  Equations  (4)  to  (6)  yields 


E 


m . V.  = M V 
11  c 


The  kinetic  energy  of  particle  i,  T^,  and  of  the  whole 
system,  T,  are  given  by 

■■’'We  use  the  following  convention:  if  a referenced  equation  appears 

in  the  same  section  as  the  reference,  then  the  section  number  is 
deleted  from  the  equation  number. 


I 


r 


_ 1 -*2  ^ 
= 2 = ? Pi  Pi 


(2-8) 


= Zt. 


(2-9) 


where'v!^  = “v.  .'v.  and^'^  =~p.  .“p. . From  Equations  (8)  and  (5) 
1 11  i 1 ^ 

we  note  that 


a T. 


a T. 


m.v.  = p. 
11  1 


1 


(2-10) 


1 

= p.  - V. 

0 p.  m. 

*^i  1 


(2-11) 


It  should  be  noted  that  in  Equation  (10)  is  considered  to  be  a 
function  of  v.  whereas  in  Equation  (11)  T.  is  considered  to  be  a 
function  of  p.  • The  confusion  between  these  two  different  func- 
tions  can  be  eliminated  by  using  different  symbols  for  the  two 
kinetic  energy  functions,  but  for  simplicity  we  will  not  do  so. 
Since  each  of  the  particles  is  independent,  Equations  (8)  to  (10) 
yield 


a T 

3 T. 
1 

a V. 
1 

a V. 
1 

^ T 

a T. 
1 

a Pi 

-> 

Pi 


(2-12) 


V 

1 


(2-13) 
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Thus,  from  the  system  kinetic  energy  function,  T,  the  i^^ 
particle  momentum  and  the  i^'^  particle  velocity  can  be  obtained 
by  partial  differentiation. 

th 

Let  f^  be  the  force  on  the  i particle.  Then  Newton's  law  of 
motion  can  be  written  as 


f t 

M 


or 


p.  = f. 
1 


-»  -» 

m.v.  = f. 
11  1 


(Z-14) 


(2-15) 


or 


*4  -» 

m.  r.  = f. 
11  1 


Equation  (14)  is  a "momentum  equation"  and  Equation  (15)  is  a 
"velocity  equation.  " Since  the  mass  m^  is  a constant,  the 
momentum  and  velocity  equations  differ  only  by  this  constant 
factor.  Let  F denote  the  total  force  on  the  system: 


(2-16) 


(2-17) 


"*  ^ E 

Now  we  can  write  f.  as  the  sum  of  f.  which  is  an  external  force 
(whose  origin  is  from  outside  of  the  deformable  body)  and  f^ 
which  is  an  internal  force  (due  to  the  coupling  among  the  mass 
points  of  the  body).  We  assume  Newton's  third  law,  that  action 
is  equal  and  opposite  to  reaction,  and  conclude  that 


-»1  ■* 

f.  = 0 

1 


(2-18) 
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Hence 


E 


i 

I 

'l 


1 1 


Combining  Equations  (6  ),  (14),  and  (19)  now  yields 
P = F = F 


(2-19) 


(2-20) 


Thus,  Newton's  law  for  a particle  also  holds  for  the  composite 
system  linear  momentum. 

Taking  the  time  derivative  of  the  kinetic  energy  of  the  i 
particle  yields 


T.  = mv.  * V.  = f. 
1 111 


(2-21) 


The  time  derivative  of  the  kinetic  energy  for  the  entire  system  is 


(2-22) 


Note  that  the  time  derivative  of  the  total  kinetic  energy  depends 
on  both  internal  and  external  forces;  i.e.  the  total  kinetic 
energy  can  change  with  time  even  in  the  absence  of  external 
forces. 

We  have  now  written  all  the  fundamental  equations  rela- 
tive to  an  inertial  origin.  In  the  next  subsection  we  write  these 
equations  relative  to  a moving  (possibly  accelerating)  origin. 


( 


I 


I 
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B.  EQUATIONS  RELATIVE  TO  A MOVING  ORIGIN 

We  now  introduce  an  arbitrary  point,  a,  with  position 

vector  r and  with  velocity  v = r . Let  R.  be  the  position 
a a a la 

vector  to  particle  i (or  point  i)  from  point  a.  From  Figure  1,  it 
is  evident  that  we  can  now  write 


r.  = r + R. 

1 a la 


(2-23) 


Multiplying  this  by  and  summing  overall  i yields 


- V 

[r  = / ->  m.  r.  = 


.r.  = Mr  + MR  (or  r = r + R ) (2-24) 

IX  a ca  c a ca' 


where 


= ^ m.  R. 


(2-25) 


Thus,  R is  the  position  vector  to  the  center  of  mass,  c,  from 

C 3i 

the  point  a.  Taking  the  inertial  time  derivative  of  Equation  (23) 
yields 


V.  = V + R. 

1 a la 


(2-26) 


The  arbitrary  point  a may  coincide  with  (be  equal  to)  c.  In  that 
case.  Equations  (23)  and  (26)  become 


r.  = r + R. 

1 C 1C 


(2-27) 
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V.  = V + R. 

(2-28) 

1 C 1C 

these  two  equations  it  follows  that 

S m.  R.  =0 

. 1 1C 

1 

(2-29) 

m.  R.  = ^ 

(2-30) 

1 1C 


Equations  (29)  and  (30)  are  consequences  of  the  fact  that  c is 
the  total  center  of  mass. 

Using  Equation  (26),  the  system  linear  momentum  becomes 


P = Mv  + MR 

a ca 


(2-31) 


where,  from  Equation  (25) 

MR  = S 


m.  R. 
ca  . 1 la 

1 


(2-32) 


If  a coincides  with  c,  then  Equation  (31)  reduces  to  P 

because  R =0. 
cc 

The  system  kinetic  energy  can  now  be  written  as 


1 -*  2 -*  ^ 1 n 2 

T = •=■  Mv  ^ + Mv  * R + y m.  R. 

2 a a ca  2 ^ 1 la 


Note  that  the  partial  derivative  of  T with  respect  to  v is  P: 


(2-33) 
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(2-34) 


3 T 


= Mv  + M R = P 
a ca 


If  a coincides  with  c,  then  Equations  (33)  and  (34)  reduce  to 


1 2 1 n 2 

T = Mv  +-5-  m.R. 

2 c Z . 


1 1C 


(2-35) 


(2-36) 


As  is  usual  when  taking  partial  derivatives.  Equations  (34)  and 
(36)  have  to  be  interpreted  carefully  because  the  two  equations 
deal  with  different  functions.  If  T denotes  the  kinetic  energy 

d ^ jJ: 

function  in  terms  of  the  independent  variable  v^,  and 
denotes  this  function  in  terms  of  v^,  then  we  can  write 

3T  * * 

O V ^ V 

a c 


From  Equations  (22)  and  (26), we  get  for  the  time  deriva- 
tive of  the  system  kinetic  energy 


where  F = F^  is  the  total  external  force.  ^ 

The  angular  momentum  about  the  point  a of  particle  i,  h^^  , 

and  of  the  whole  system,  H , are  given  by 

a 
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h . = R.  X p. 
ai  la  1 


il  =E 


a . ax 
1 


(2-39) 


Thus,  H is  given  by 

3. 


H = 2 R;  X P:  = S 


-»  -> 

m.  R.  X V. 
1 la  1 


(2-40) 


K a coincides  with  c,  this  expression  becomes 


H„  = S R.^  ^ P:  = ^ 


I^=  2^  Ri^x  Pi=  ^ "^i^^ic”'  '^i 
i 1 


(2-41) 


If  we  write  R.  = R.  + R in  Equation  (40),  we  find  that 
la  1C  ca 


H = H + R X P 
a c ca 


(2-42) 


We  get  an  alternate  expression  for  H by  expressing  in  terms 


of  V as  given  by  Equation  (26).  Then 
3 

-♦  ^ 

H = MR  XV  + m R X R 

a ca  a x la  la 


(2-43) 


If  a coincides  with  c,  we  get  an  alternate  expression  for 


H = / m.  R.  X R. 
c . X xc  xc 
X 


(2-44) 


Taking  the  time  derivative  of  Equation  (40)  yields 
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I 


-> 

H 

a 


2 R.  X p.  + 
. la  1 


(2-45) 


We  substitute  p.  = f.  on  the  right  hand  side,  and  then  introduce 
the  torque  or  moment  of  force,  L^,  on  the  system  about  the 
point  a: 

->  -t  ^ -»  -» E -»  E 

L = R.  X f.  = Z-»  R.  X f.  = L ^ (2-46) 

a . la  1 . la  1 a 

1 1 


where  we  have  assumed  that  action  is  equal  and  opposite  to 
reaction  and  hence 


Z 


R.  X f/  = 0 
la  1 


(2-47) 


Next  we  note  that 


^ia  Pi  = V (^i  ■ ^a>  * Pi  = -'^a  " ^ 
i 1 

Combining  Equations  (45)  and  (48)  and  making  use  of  Equation  (47) 
now  yields 


H + V X P = L 
a a a 


(2-49) 


K a coincides  with  c,  this  reduces  to 


since  v x p = 6 , 
c 


(2-50) 


4 
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'-'"HI"  "■  ill. 


The  definition  of  the  total  torque  about  the  point  a. 

Equation  (46),  and  the  definition  of  the  total  angular  momentum 
about  the  point  a.  Equation  (40),  are  very  similiar.  This  simi- 
liarity  can  be  enhanced  by  introducing  the  moment  of  force  about 
the  point  a on  particle  i,  as  follows 

t . = R.  X ?.  (2-51) 

ai  la  1 

Then 


Note  the  similiarity  of  Equations  (51)  and  (52)  to  Equations  (38) 
and  (39),  respectively. 

The  angular  momentum  is  a fundamental  quantity  but  we 
cannot  obtain  it  by  partial  differentiation  of  the  kinetic  energy 
function  given  in  Equation  (33).  This  situation  will  be  remedied 
in  the  next  subsection  where  we  introduce  an  angular  velocity,  and 
then  the  angular  momentum  can  be  obtained  by  partial  differen- 
tiation of  the  kinetic  energy  with  respect  to  this  angular  velocity. 

C.  EQUATIONS  RELATIVE  TO  A ROTATING  FRAME 

In  the  last  subsection,  we  obtained  the  following  expressions 
for  the  linear  momentum,  P,  angular  momentum  about  a,  H^, 
kinetic  energy,  T,  and  time  derivative  of  kinetic  energy,  T: 

P = Mv  + ^ m.  R.  (2-53) 

a . 1 la 

1 


% 

»: 


4 
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i 


= MR  ; 

-> 

< V + 

Z!  m R. 

•-* 

X R. 

(2- 

■54) 

ca 

a 

1 la 
1 

la 

1 -»  2 
: ^ Mv 

4 -> 

m.  R.  * V 

4 

m.  R. 

(2- 

■55) 

2 a 

i 

1 la  a 

2 

. 1 la 

1 

-»  -» 

: F • V 

a 

- 2: 

f.  • R. 

1 la 

(2- 

-56) 

We  will  now  replace  the  inertial  time  derivative  of  in 

these  equations  by  the  time  derivative  with  respect  to  a frame 

-♦B 

(which  we  can  call  "frame  B")  which  has  angular  velocity  ^ 
with  respect  to  inertial  space.  We  know  that  for  any  vector  V 
we  have 


B, 


V = V + u)  X V 


(2-57) 


where  V denotes  the  time  derivative  of  V with  respect  to  the 
frame  with  angular  velocity  u)  (i.  e.  with  respect  to  frame  B). 
Applied  to  the  vector  Equation  (57)  yields 

i*  S ->B  -> 

R.  = R.  + u)  * R. 
la  la  la 


(2-58) 


-*B  ^ 

It  is  convenient  to  let  u.  denote  R. 

la  la 

-♦B  - I 
u . = R. 

la  la 


(2-59) 


It  is  also  convenient  to  write 

*B 


-»R  -»  ~t 

X R.  = R 

la  la 


U) 


(2-60) 
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where  R. 

vector  R.  , 
la 

R.  . Thus, 
la 


is  the  skew- symmetric  dyadic  formed  from  the 

and  R^  is  the  dyadic  transpose  (or  conjugate)  of 

R.  can  be  expressed  as 
la 


, ST 
u.  + R. 
la  la 


(2-61) 


Substituting  Equation  (61)  into  Equations  (53)  to  (56)  yields 


~t  -*B 

P=Mv  + MR  • UJ  + Z-.  m u 

a ca  ^ 1 la 


(2-62) 


-»  . . :^B 


a 

— AVI 

ca 

a a 

1 

1 

la 

T = 

1 -»  2 
^ Mv 
2 a 

V -►B 

+ m.u  . 

1 la 

1 

-> 

• V + V 

a a 

mr' 

4 

-*B  ^ 

CD  • 1 

a 

• (i)  + Z-rf  m 

• R^ 

1 la  la 

*B 


t . 

ca 


(2-63) 


(2-64) 


T = F * V + L 
a a 


-►B  V*  ■?  ■*£ 

• f.  • u . 


-+B 
u . 

1 la 


(2-65) 


where 


4^ 

I 

a 


m.R. 


1 la 


(2-66) 


Evidently,  is  the  inertia®^  of  the  deformable  body  about  the 
arbitrary  point  a.  If  a coincides  with  c,  we  have 


See  Appendix  A 
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^ ,, 


7 = E 


m.R.  • R !■ 

1 1C  1C 


(2-67) 


If  we  write  R.  — R.  R in  Eciuation  (66),  we  find  tliat 
la  1C  ca 


I + MR  • R 
a c ca  ca 


(2-68) 


Considering  v and  u5  ® to  be  independent  variables  in  the 
a 

kinetic  energy  function  given  by  Equation  (64),  we  see  that 


3 T 


= 


(2-69) 


3 T 


B = «a 


(2-70) 


Here  the  partial  derivatives  of  T are  understood  to  mean  the 


partial  derivatives  of  the  function  T given  in  Equation  (64). 

The  time  derivatives  of  P and  H are  still  given  by 

a 


P = F 


(2-71) 


H + V 
a a 


P = L 


(2-72) 


Equations  (69)  to  (72)  are  very  simple  but  very  general  equa- 
tions which  hold  for  any  deformable  body,  even  a body  containing 
fluids®^.  However,  Equations  (69)  to  (72)  describe  only  the  6 
external  degrees  of  freedom  of  the  deformable  body.  In  order  to 
describe  the  n ^ 3N  - 6 internal  degrees  of  freedom  of  the  system 


-23- 


of  N particles,  we  must  somehow  introduce  n appropriate 
coordinates  and  velocities.  In  addition,  it  would  be  desirable 
to  introduce  n appropriate  momenta. 

It  is  instructive  to  write  Equations  (62)  and  (63)  in  the 
following  matrix  form 


p 

mS 

MR^ 

ca 

4* 

m^^E 

4^ 

m^E 

44  n 

H 

a_ 

MR 

ca 

t 

a 

m R 
1 la 

2 2a 

”"N^Na. 

♦B 

^la 

^2a 


Na 


where  E is  the  identity  dyadic.  Since  the  system  has  only  n 

internal  degrees  of  freedom,  the  relative  velocity  vectors 

-*B  -»B  ~*B 

u 1^,  u • * ■ , u must  be  expressed  in  terms  of  only  n 
coordinates  and  velocities.  Also,  there  are  n more  momenta 
which  characterize  the  internal  degrees  of  freedom  of  the 
system. 


See  Appendix  B. 
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m.  INTERNAL  GENERALIZED  COORDINATES 


If  there  are  n internal  degrees  of  freedom  described  by 

internal  generalized  coordinates  ^ an’ 

u?  must  be  a linear  combination  of  the  time  derivatives  of  these 
la 

internal  generalized  coordinates.  Thus,  we  can  write  an  equation 
of  the  form 


-►B 

u . = 


la 


i = 1 


(3-1) 


where  . in  general,  depends  on  the  internal  generalized 
B B B 

coordinates  £ i»  C ■>«"'*>  K however,  ?.  does  not  depend 

on  the  time  derivatives  of  the  internal  generalized  coordinates. 

Now  R.  is  a function  of  the  internal  generalized  coordinates, 
la  ^ 

and  we  assume  that  R.  is  not  an  explicit  function  of  time.  Hence 

la 


- R 
u . - R. 

la  la 


3R.  .X. 

13  ^ B 


aj 


(3-2) 


Comparing  Equations  (1)  and  (2),  we  see  that 

3r. 

la 

a 

aj 


(3-3) 


Equation  (1)  can  be  written  conveniently  in  matrix  form  as 
follows  ; 


See  Appendix  B 


(3-8) 


Equation  (1)  can  now  be  written  as 


-►B  xB^  :b  ;b^  xB 

u . = * . £ = £ t . 

la  la  ’ a a la 


(3-9) 


**B 

These  two  different  forms  of  expressing  u.  are  useful  in  form- 
^b2 

ing  quantities  like  u^^  : 


-»B^  ->B  . -*B  ;B^  xB  xB^  ;B 

u.  =u.  u.  =£  * £ 

la  la  la  a la  la  a 


(3-10) 


where  . ^B 

is  an  n X n matrix  of  scalars: 


la  la 


. $ 


B 


la  la 


f?B  |B  ...  1 

• L ial  iaZ  ianj 


(3-11) 


lan 


tB 

tB  . 

lal 

lal 

lal 

iaZ 

ial 

ian 

.pi  ... 

, • 

laZ 

lal 

laZ 

iaZ 

laZ 

lan 

lan 

lal 

lan 

. ... 

iaZ 

■ 

lan 

lan 
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Thus,  the  element  in  the  row  and  column  of  ‘ is 


V la  xa 


=tB  . 

a/jk  laj 


(3-12) 


Equations  (9)  and  ( 1 0 )can  now  be  used  to  write  Equations  (2-62) 
to  (2-65)  in  the  following  form 


P = Mv  + MR 

a ca 


1 -*  2 •* 
i-  Mv  + V 
2 a a 


-♦B  *B^  ;b 

1 la  a 

(3-13) 

1 

+ y • 3B+  Zm.R.  • 

a . 1 la  la  ^a 

1 

(3-14) 

~t  "*B  y^  B^  iB 

• MR  • uu  + V • ^ m4  i 

ca.  a . 1 la  a 

(3-13) 

1 -♦B  ♦♦  -►B , -»B 

_ (1)^  - I *00  + u)  • m. 


■ R.  • 

1 la  la  a 


4 p' 

1 


Z ?.  ■ 


T = F*v+L*uJ  +il-<f.''*.  ^ 

a a ^ 1 la  a 


(3-16) 


Note  that  in  each  of  these  equations  we  have  a summation  over 
all  the  particles  (from  i = 1 to  N).  We  can  eliminate  this 
explicit  dependence  on  the  individual  particles  by  introducing 

’’A 

the  following  notation 


= Z 


m.  s . 

1 la 


(3-17) 


* See  Appendix  B 
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1 


These  equations  can  be  written  conveniently  in  matrix  form  as 
follows  : 


H MR  I 

a ca  a 


T 4 e 

2 a ^ a 


ME  MR 


ca  a 


T = F L 

a a _ 


In  addition,  we  have  the  relationships 


* See  Appendix  B 


(3-28) 


(3-29) 


(3-30) 


(3-31) 


I 


(3-32) 


St  •+ 

0 0) 


where  here  the  partial  derivatives  of  T mean  the  partial  deriva- 

^ * 3 

tives  of  the  function  of  v , , and  £ (and  of  the  internal 

^ B B ^ B 

generalized  coordinates  £ i » ^ 4 ) given  in  Equation 

di^  3iTi 

(29). 

We  have  now  achieved  our  objective  of  expressing  the 
fundamental  quantities  in  terms  of  the  velocities  of  the  6 external 
degrees  of  freedom  plus  the  velocities  of  the  n internal  degrees 
of  freedom.  However,  so  far,  we  only  have  the  external  momen- 
ta Pand  i?  ; we  would  also  like  to  have  n internal  generalized 
s, 

momenta.  Since  the  external  momenta  can  be  obtained  by 
partial  differentiation  of  the  kinetic  energy  with  respect  to  the 
angular  velocities,  we  are  led  to  define  the  n internal  generalized 

• B 

momenta  by  partial  differentiation  of  T with  respect  to  £ . De- 
noting  the  internal  generalized  momentum  by  g we  thus  have 


B _ 


X rB 
V + P 

a a 


a ^ a 


(3-33) 


Combining  Equations  (28)  and  (33)  now  yields 


(3-34) 
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T = — 

2 


♦B 


II 

P H 

a 

^a  J 

• 

p 

1 

H 

a 

B 

Sa 

— — 

(3-35) 


We  now  have  a full  set  of  external  and  internal  velocities, 

plus  a full  set  of  external  and  internal  momenta.  A full  set  of 

^ ^ . B 

dynamics  equations  is  a set  of  equations  for  P,  H , and  g , or 

s*  .‘g  ^ ^ 

a set  of  equations  for  v , , and  ^ . Equations  (2-71)  and  (2-72) 

a^»  a 

give  us  expressions  for  P and  H , respectively:  this  will  be  a 

full  set  of  dynamics  equations  as  soon  as  we  get  an  equation  for 

g (see  Section  V,  page  65). 

^ ^ ^ -> 

Recall  that  we  defined  P and  in  terms  of  p^  (Equations  (2-6) 

and  (2-40]^espectively).  Thus,  we  were  able  to  obtain  equations 

■4  4 

for  P and  H by  simply  differentiating  these  defining  equations 
^ 4 

with  respect  to  time,  and  then  using  Newton's  law  to  replace  p. 

-♦  B ^ 

with  f . . In  the  next  section  we  will  define  g as  being  given  by 


ia 


-»  . B 

p..  From  this  definition  we  will  obtain  g as 


Z 


$ ? • p.  + S '5’® 

la  la 

1 1 


p. . When  we  replace  p.  with  f.  we 

B 


notice  that  the  second  of  these  summations  is  simply  k as 

a 

defined  in  Equation  (22).  We  will  then  have  our  equation  for 

- B . V jB  . ^ , 

e as  soon  as  we  put  / . « . p.  in  a more  convenient  form 
^ a la  '^i 

1 

which  does  not  involve  an  explicit  summation  over  particles.  Of 

course,  we  also  have  to  show  that  the  definition  of  g^  as 

xB  B ^ 

/ x • p.  is  equivalent  to  the  definition  of  g as  in 

i 


Equation  (33). 
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i 


IV.  TRANSFORMATION  OPERATOR  EQUATIONS"' 


We  will  now  use  the  transformation  operator  approach  to 
derive  the  desired  full  set  of  equations  of  motion  for  our  deform- 
able body.  The  transformation  operator  formalism  has  been 
documented  in  References  (66)  and  (68)  and  therefore  we  will  give 
here  only  a brief  summary  in  a form  most  suitable  for  applica- 
tion to  a deformable  body.  It  is  interesting  to  note  that  the  en- 
tire formalism,  which  is  based  on  Kron's  method  of  subspaces 

is  essentially  a generalization  of  linear  matrix  structural  analy- 
83-84 

sis.  The  generalization  consists  mainly  of  starting  out 

with  coordinate  dependent  velocity  transformations  rather  than 
with  linear  and  constant  coordinate  transformations.  In  general, 
coordinate  transformations  are  nonlinear,  but  the  corresponding 
velocity  transformations  are  always  linear. 


A.  BASIC  FORMALISM 

The  equations  of  motion  for  a dynamical  system  with  a 
finite  number  of  degrees  of  freedom  can  be  put  in  the  form 


G + X = K (4-1) 

where  G is  the  system  momentum,  K is  the  system  force,  and 
X is  an  extra  term  which  is  quadratic  in  the  momentum.  In 
general  X also  depends  on  the  coordinates.  The  momentum  G is 
linearly  related  to  the  system  velocity  o : 

G = y.  • cf  (4-2) 


where  M-  is  the  system  mass.  Taking  the  time  derivative  of 

* The  final  equations  obtained  in  this  section  are  summarized  in 
Section  V,  page  65.  The  reader  may  wish  to  look  at  Section  V 
before  reading  Section  IV. 
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Equation  ( 2 ) and  then  substituting  into  Equation  (1)  yields 


'I 


M.  • a + Y = K 


(4-3) 


where 


Y = p,  * a + X (4-4) 

Equation  (1)  is  a "momentum  formulation"  whereas  Equation  (3) 
is  a "velocity  formulation.  " The  principal  difference  between 
the  momentum  and  velocity  formulations  is  that  in  the  momentum 
formulation  we  have  the  extra  term  X whereas  in  the  velocity 
formulation  we  have  the  extra  term  Y.  Often  it  is  simpler  to 
obtain  X than  it  is  to  obtain  Y (this  is  certainly  the  case  if  to 
obtain  Y,  we  first  obtain  X as  indicated  in  Equation  (4) 

In  addition  to  Equations  (1)  to  (4),  we  also  have  the 
kinetic  energy  expressions 


T = 


(4-5) 


• t 

T = K • a 


(4-6) 


5 T 

3 CT 


= G 


(4-7) 


When  performing  the  partial  differentiation  of  T with  respect  to 
c we  consider  T to  be  the  quadratic  function  of  a given  in 
Equation  (5). 

The  system  mass  is  positive  definite  and  symmetric,  and 
therefore  it  has  a positive  definite  symmetric  inverse  : 


1 


1 

1 

1 


- H- 


V 


(4-8) 


-1 

p,  = 


Therefore,  Equation  ( 2 ) can  be  inverted  as  follows 


a = V • G 


(4-9) 


Consequently,  the  kinetic  energy  can  also  be  expressed  as 


T = 


V • G 


(4-10) 


From  this  follows  the  relationship 
3 T 

-2-  = CT  (4-11) 

3G 

where  here  we  consider  T to  be  the  quadratic  function  of  G given 
in  Equation  (10), 

We  now  make  a linear  velocity  transformation 

a = A • CT  (4-12) 


i 

1 

I 


where  o is  a new  system  velocity  and  A is  the  transformation 
operator.  In  conjunction  with  this  velocity  transformation,  we 
also  make  the  following  momentum  and  force  transformations 


G = a‘  • G 
K = A*  * K 
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(4-13) 

(4-14) 


„ 


We  also  make  a congruence  transformation  on  the  mass: 


M,  = • M«  * A 

Substituting  Equation  (12)  into  Equations  ( 5 ) and  (6)  and  then 
making  use  of  Equations  (13)to  (15)yields 

T =La^  • m:  • a =jC^  ' a 
t = • CT 

Thus  CT,  Tf,  G,  and  K are  properly  defined  in  order  to  keep  the 
functional  form  of  T and  T invariant. 

Differentiating  G with  respect  to  time  yields 

G + X = K (4-18) 

where 

X = A^  • X - ’ G (4-19) 

We  also  find  that  G and  c are  related  as  follows 
— _ _ S T 

G = ij,  • a = (4-20) 

9 a 

Substituting  Equation  ( 2 0 )into  Equation  (18)  yields 

• ? + Y = K (4-21) 


where 


Y = M.  • a + X 


(4-22) 


We  get  an  alternate  expression  for  Y by  substituting  Equation  (12) 
into  (3)  and  then  multiplying  the  resulting  Equation  (3)  from 
the  left  by  a'".  The  resulting  Y is 

Y=A^‘Y  + A**M.‘A’CT  (4-23) 

If  Equation  (22)  is  used  to  obtain  Y,  then  it  obviously  requires 
more  effort  to  obtain  Y than  to  obtain  only  X.  However,  if 
Equation  (23)  is  used  for  Y,  it  is  not  clear  whether  it  is  simpler 
to  find  this  Y or  the  X given  by  Equation  (19).  In  fact,  by 
writing  G = M>*a=M'*A  *awe  can  put  X in  the  form 

X = A'^*X-A’^-p-A-a  (4-24) 

Comparing  Equations  (23)  and  (24),  we  note  that  both  X and  Y 

appear  to  be  equally  complicated. 

In  general,  the  transformation  operator  A is  not  invertible, 

but  it  is  always  one-to-one  (alternatively,  it  always  has  full  column 

8 5 *■  8 6) 

rank).  Therefore,  A has  a left  inverse  B so  that 


B • A = 1 (4-25) 

where  1 is  the  identity  of  the  same  dimension  as  a (and  also  of 
the  same  dimension  as  G,  K,  X,  and  Y).  From  Equation  (25) 
it  follows  that 


I 


I 
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(4-Z6) 


i 

j 

i 


Thus,  A*  has  a right  inverse  However,  in  general,  A does 

not  have  a right  inverse,  and  therefore  A*'  in  general  does  not 
have  a left  inverse.  Consequently,  Equation  (12)  can  be  solved 
for  a : 


a = B • a (4-27) 

but  Equations  (13)  and  (14)  cannot  similiarly  be  solved  for  G 
and  K in  terms  of  G and  K,  respectively.  Similiarly,  Equation 
(15)  cannot  be  solved  for  g.  in  terms  of  U. 

Since  y.  is  positive  definite  symmetric  and  since  A is  one-to- 
one,  it  follows  from  Equation  (15)  that  TI  is  also  positive  definite 
symmetric.  Therefore  F has  a positive  definite  symmetric  in- 
verse V : 


_-l 


= V 


(4-28) 


Thus  <T  can  be  expressed  as  follows: 

a T 

a = V . G = a G 


It  can  be  shown  from  general  tensorial  considerations  that 
there  exists  a C such  that 


■q;  = C • TI  + H • C*^  (4-30) 

X = -C  • G = -C  • p;  • c (4-31) 


* C is  linear  in  a and  in  the  Christoffel  symbols  (which  general- 
ly are  functions  of  the  coordinates). 


I 


I 

-.1 
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(4-32) 


Y = M.  • • a 

Now  define  D by 

D = C * H (4-33) 

Then  Equations  (30)  to  (32)  become 


Tl  = D + 

(4-34) 

XI 

II 

Dl 

Q1 

(4-35) 

Y = • o 

(4-36) 

where  we  have  made  use  of  the  fact  that  p.  is  symmetric 

(T  = il*').  From  these  equations  it  is  particularly  evident  that 

Y - X = fr  • a (4-37) 

From  Equations  (35)  and  (36)  it  is  evident  that  if  we  obtain  X 
and  Y by  first  obtaining  D,  then  X and  Y are  obtained  with  equal 
difficulty. 

B.  APPLICATION  TO  A DEFORMABLE  BODY 

In  order  to  use  the  transformation  operator  formalism,  we 
must  start  out  with  a set  of  O',  g.,  G,  K,  X,  and  Y which  satisfy 
Equations  (1)  to  (3).  Then  we  must  introduce  a new  velocity 
'o  and  a transformation  operator  A.  The  transformation  opera- 
tor formalism  then  does  the  rest. 
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We  consider  a system  of  N particles  and  we  let  a be  a 

th  "* 

column  matrix  of  N vectors,  whose  i element  (vector)  is  v^, 
the  velocity  of  the  i^^  particle: 


'N 


(4-38) 


We  define  p to  be  an  N x N matrix  of  dyadics  which  is  diagonal 

th  ** 

and  whose  i diagonal  element  is  m^E  , where  m^  is  the  mass  of 
the  i*^^  particle: 


M-  = 


mj6 


m_E 


(4-39) 


where  the  off-diagonal  elements  of  p.  are  the  zero  dyadic.  Since 
G = p ' cT  we  see  that  G must  be  a column  matrix  of  N vectors, 
whose  i*'^  element  is  p^,  the  linear  momentum  of  the  i*’^  particle: 
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G = 


(4-40) 


Pn 


Next,  we  define  K to  be  a column  matrix  of  N vectors,  whose 
i^'^  element  is  L , the  force  on  the  i''^  particle: 


(4-41) 


Since  we  have  = L we  obviously  have  X equal  to  zero: 


X = 0 (4-42) 

Thus,  X is  a column  matrix  of  N zero  vectors.  Similiarly, 

since  m.  v.  = f.  , we  also  have  Y equal  to  zero: 

11  1 

Y = 0 (4-43) 
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i 


g 


and  therefore  Y is  also  a column  matrix  of  N zero  vectors. 

Thus,  the  equations  of  motion  for  the  system  of  particles  are 


G = K 


(since  X = 0) 


(4-44) 


G = li.  • a 


(4-45) 


These  equations  represent  the  momentum  formulation.  The  velo- 
city formulation  equation  is 


M.  • a = K 


(since  Y = 0) 


(4-46) 


We  now  define  a as  follows; 


CT  = 


-> 

V 


a 

:b 


(4-47) 


where  v is  the  inertial  linear  velocity  of  an  arbitrary  point  a as 
introduced  in  Section  II  , id  is  the  inertial  angular  velocity  of  an 


:B  . 


arbitrary  frame  B as  introduced  in  Section  U , and  4 is  an  n 

th  * B ^ 

element  column  matrix  whose  j element  is  4 • » which  is  the 

th 

time  derivative  of  the  j internal  generalized  coordinate. 


We  determine  the  transformation  operator  A by  combining 
the  following  three  equations  (see  Equations  (Z-26),  (Z-6l),  (3-1),  and 
(3-9)) 


v.  = v + R. 

1 a la 


(4-48) 


-4Z- 


Ji 


-♦  -*B  t -♦  B 

R.  = u + R.  • 
la  xa  la 


-B  ^ 

la  laj  ^aj 


la  a 


Thus  we  have 


"*  "*  B , iB  »B 

V.  = v+R.  ’UJ  +'^■4, 
1 a la  la  a 


= E R.^ 

la  la 


Thus,  A is  given  by 


t* 

~t 

B^ 

E 

Bi 

la 

'la 

~t 

B^ 

E 

^2a 

' 2a 

E R^ 

^ ^Na  Na 


G is  now  defined  by 


L 


G = A • G 


4* 

E 


4» 

E 


la 

2a 

^Na 

^a 

‘L  - 

Na 

Pi 

P2 


(4-53) 


N 


E 


Pi 

R.  • 

P- 

la 

-4 

P- 

la 

From  Equation(2-6)we  see  that  the  first  element  of  G is  P,  the 

total  system  linear  momentum.  From  Equation(2-40)we  see  the 



second  element  of  G is  H , the  total  system  angular  momentum 

^ — B 

about  a.  We  will  denote  the  third  element  of  G by  g . Later,  we 

3l 

will  see  that  this  is  the  same  quantity  which  we  introduced  in 
Equation{3 -33)  and  which  we  called  the  internal  generalized  momen- 
tum. Thus,  G is  given  by 


(4-54) 
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K is  obtained  similarly; 


Zt 


K = A • K 


R.  • f. 
la  1 


(4-55) 


where  we  have  made  use  of  Equations( 2-17)^2 -46),  and(3-22).  Thus, 
the  first  element  of  K is  the  total  (external)  force  on  the  system, 
the  second  element  of  K is  the  total  (external)  torque  on  the 
system  about  a,  and  k is  a quantity  which  we  can  call  the  force 

on  the  internal  degrees  of  freedom.  If  we  split  f.  up  into  an  ex- 

-»E  ^ 

ternal  part  f.  (whose  origin  is  from  outside  of  the  deformable 

^ ■*! 

body)  and  an  internal  part  f^  (duf>  to  coupling  among  the  mass 
points  of  the  body),  then  we  can  write 


, B , BE  , , BI 
k = k + k 
a a a 


(4-56) 


where 


la  1 


= E • ? 


(4-57) 


(4-58) 


P is  obtained  as  follows 
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n'l . E m . R . 

1 . 1 la 

1 1 


Z— < m . s . 


M.  = A • n ' A = 


m.R.  i—t  m.R.  R. 

1 la  . 1 la  la 

1 1 


m.R.  • $ 

1 la  la 


m . $ m . $ ^ • R 


1 la  . 1 la  la 

1 1 


V .B  . xB* 
/ -»  m.  ? . • ? . 


1 la  la 


where  we  have  introduced  M from  Equation(2-3)^MR  from 

Equation(2-25)jjI  from  Equation(2-66),a  from  Equation(  3-17), 

B ^ B ^ 

0 from  Equation(  3 -19),  and  y from  Equation(3 -21) . 
a a 

The  equation  0=11  " ^ is  now  precisely  Equation(  3-34) 

given  in  Section  III.  This  proves  the  statement  that  g as  defined 

a 

in  this  section  is  the  same  quantity  as  defined  in  Section  III: 


If  the  arbitrary  point  a coincides  with  the  total  center  of  mass  c, 

then  P=  Mv  and  R =0  imply  that 
c cc  ^ ^ 


.Bx  E 


m . i ? =0 

1 1C 


L 


Hence,  in  this  case,  the  internal  generalized  momentum  is 
given  by 


B 

8c 


Pi 


,B 


+ y® 
c 


(4-62) 


where 


3® 

c 


E 


m.  9. 

1 1C 


5 B 

m.  9 . 

1 1C 


ic 


(4-63) 


(4-64) 


Note  that  after  the  arbitrary  point  a has  been  specified  (say,  by 

• B 

letting  a coincide  with  c),  we  still  have  arbitrariness  in  g^ 

and  £ due  to  the  arbitrariness  of  the  frame  B. 
a 

C.  THE  FORM  OF  X AND  Y FOR  A DEFORMABLE  BODY 

So  far,  our  application  of  the  transformation  operator  form- 
alism to  a deformable  body  has  not  really  yielded  much  more 
than  could  be  obtained  from  linear  matrix  structural  analysis. 
However,  it  should  be  noted  that  our  formalism  has  additional 
flexibility  because  we  are  actually  using  a coordinate  dependent 
transformation  matrix  A [Equation  (5Z)  ],  and  this  transforma- 
tion matrix  has  some  elements  which  are  dyadics  and  some 
elements  which  are  row  matrices  of  vectors.  The  key  ingredients 
in  a set  of  exact  dynamics  equations  is  the  extra  term  X or  Y, 
which  depends  quadratically  on  the  momentum  or  velocity.  The  X 
and  Y terms  cannot  be  obtained  from  linear  matrix  structural 
analysis  because  these  terms  depend  on  the  time  derivative  of 
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X = -A*  • G = -A*  • M.  • A • a 


(since  X = 0)  (4-65) 


Y = A • ti.  • A • cj 


(since  Y = 0)  (4-66) 


Comparing  Equations  (65)  and  (66)  with  Equations  (34)  to  (36), 
we  see  that  in  this  case 


D = A • \i  • A 


(4-67) 


S m.R. 


m.R.  • r: 

1 la  la 


m.R.  • 

1 la  la 


S m.i?  ^ m.l?  • R^  X]  m.l?  • 


1 la  la 


1 la  la 


where  O is  the  zero  dyadic,  and  0 represents  a row  matrix  of  n 
zero  vectors.  We  can  write  D more  compactly  as  follows: 


(4-68) 


the  transformation  operator  A.  We  now  turn  to  determining  the 
form  of  X and  Y for  the  deformable  body. 

Since  in  the  present  application  both  X and  Y are  zero,  we 
have  from  Equations  (23)  and  (24) 
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whe  re 


i.R.  • R. 

1 la  la 


• R.^ 


1 la  la 


= ^ m.R.  • 


1 la  la 


1 la  la 


(4-69) 

(4-70) 

(4-71) 

(4-72) 


°22  ^22 


°32  °23 


°23  °32 


^33  ■*'  ^33 


Taking  the  time  derivative  of  >1  yields 


** 

~t 

• B 

o 

MR 

a 

ca 

a 

0 

ab' 

MR 

I 

3 

ca 

a 

a 

yf 

a 

a 

+ we  see 

that  we  now 

have 

t 

ft 

°22  ^ ^22 

= I 

a 

(4-75) 
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= D„  + D,  * 
a 3Z  Z3 


°33  ^ °33 


(4-76) 

(4-77) 


We  can  now  write  X and  Y as 


X = -D  • a = - 


Y = • a 


MR  • V + D_,  • + D_, 

ca  a ZZ  Z3  a 


a®  • V + D • + D 

a a 3Z  33  a 


ca  a a 


D ‘ + D ^ 

ZZ  3Z  a 


D ^ + D * 

Z3  33  ^a 


(4-78) 


(4-79) 


We  will  now  show  that  the  second  element  of  X is  v x P as 

a 

required  by  Equation(Z-49).  Note  that 


-MR  ' V = -M(v  - V ) X V - -Mv  X V 
ca  a c a a c a 


(4-80) 


= V X P 
a 


Next  we  show  that 


D22  • 0 


Z3  ^a 


(4-81) 
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From  Equations  (^9)  and  (71),  we  see  that  this  requires 


m.R. 

1 la 


~t  -»B  ~ 

R.  • (JO  + 2-»  r^-R- 
la  ■ 1 la 

1 


5®'  s'"  = 0 

la  a 


(4-82) 


Factoring  out  the  common  term  m.R.  and  introducing  R from 
® 1 la  ® la 

Equations  (48)  and  (51)  allows  us  to  write  Equation  (82)  as 


S m R 


1 la 


R.  = 0 
la 


(4-83) 


This  equation  is  obviously  satisfied  and  therefore  Equation  (81) 
is  established.  Thus,  Equation  (78)  can  be  written  as 


X-* 

V 


-♦B 

U) 


:B 


MR  • V 
ca  a 


a®  • V + D • 
a a 32  33  ^ a 


V X P 

a 


;B 


(4-84) 


where  we  have  introduced  the  symbols  X~*  , X.^  , and  X „ for 

''a 


€ 


the  elements  of  X.  Note  that  X „ is  an  n element  column  ma- 

th  ^ 

trix  of  scalars:  the  j element  of  this  matrix  is  denoted  by 


aj 


The  expression  for  Y given  in  Equation  (79)  cannot  be  sim- 
plified. However,  we  can  get  an  alternate  expression  for  Y from 
Y = TI  • a + X.  Using  Equation  (74)  for  II  and  Equation  (84) 
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for  X yields 


Y = 


or 


Y = 


MR  • il®  + a® 

ca  a a 

1 

to 

1 

~ -*  p,  • 

MR  ■ V + I • ir  “ + 0 

ca  a a 

b’^  aB 
a ^ a 

+ 

-MR  • V 

ca  a 

O'®  * V + 8®  • (T®  + 

a.  a a a 

1 

m 

■ UJ> 

IX 
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— 

Y-* 

V 
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0) 


Y-r 

L 


MR  ^ 

ca  a a 


*?  • ^ .B 

a a ^ a 


D ^ + D 

23  ^33  ^a 


(-1 


(4 


where  we  have  introduced  the  symbols  Y**  . Y_,„  , and  Y „ for 

a 

the  elements  of  Y.  Y „ is  an  n element  column  matrix  of  sea- 

«a 

th  ^ 

lars;  the  i element  of  this  matrix  is  denoted  by  Y ^ . 

k . 

aj 

Lagrange's  equation  for  the  internal  generalized  coordinates 
takes  the  form 


d /a  T 
dt  5^ 


St  , b 

— =^a 


(4 


-85) 


-86) 


-87) 


w 


The  first  term  on  the  left  is  Therefore,  the  second  term  on 

a _ 

is 


the  left  is  X.  . Thus,  an  alternate  expression  for  X 


X 


a T 


B 


a? 


B 


(4-8H) 


We  will  return  to  this  equation  after  we  get  explicit  expressions 
for  D and  X. 


D.  EXPLICIT  EXPRESSIONS  FOR  D,  X,  AND  Y FOR  A 


DEFORMABLE  BODY 


The  expressions  for  D,  X,  and  Y given  above,  show  the 


general  form  of  these  quantities,  but  they  do  not  show  explicitly 


the  dependence  of  these  terms  on  the  velocities  v , and  4^  . 

^ 3. 


From  Equation  (67),  we  note  that  in  order  to  get  explicit  exprcs- 

— ~ • B "*B 

sions  for  D,  we  need  to  express  R.  and  $ . in  terms  of  v , a , 


and  4 


:B 


la 


la 


From  Equations(  Z-58),  p-59^(  3-1),  and(3-9)wc  have 


I =4® 

la  a aa 


R. 

la 


(4-89) 


and  from  this  follows 
,t 


R.  =4"" 
la  a xa 


R. 

la 


R. 


-B 


(4-90) 


^ B B 

where  ^ . is  formed  from  $ . by  replacing  all  vector  elements 
la  la  ^ 


by  the  skew- symmetric  dyadic  of  these  vectors.  Thus, 


1 
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ial 
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'I'® 

' iaZ 


(4-91) 


ian 


Taking  the  transpose  of  Equation  (90)  yields 
,t 


R.^  = - R.^ 

la  la  a la 


ujB  ■^B 


R. 


t 


la 


(4-92) 


where  we  have  changed  signs  in  the  second  and  third  terms  by 
using  a)  = -u)  . It  should  be  noted  that  $ is  the  transpose  of 
of  Equation  (91): 


la 


la 


r t 

t 

'I'B 

lal 

laZ 

^ ial 

laZ 

t 

^ . we 

first 

la 

matrix, 

and  V 

lan  J 


3 

ian 


(4-93) 


rB 

ia 


rBt 


replace  each  clement  by  its  transpose.  We  do  not  form  by 

B^ 

first  forming  i . and  then  replacing  each  vector  elements  by  its 
la 

symmetric  dyadic;  in  fact,  as  shown  in  Equation  (93),  we 
actually  want  the  negative  of  this.  It  is  easy  to  see  that 
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(4-94) 


^ a la 


la  a 


Comparing  Equations  (90),  (92),  and  (94)  now  shows  that  this 

peculiar  minus  sign  is  just  what  is  required  to  get 

Since  is  a column  matrix  of  vectors,  we  get,  by  a simple 
la 

generalization  of  Equation(  2-57),the  result 


a ia  la 


B 

where  is  the  time  derivative  of  with  respect  to  frame  B: 

la  la 


® = E 

1 n * a 


ak  3 .B 
^ ak 


It  is  now  convenient  to  introduce  the  operator  „ as  follows 


^^al 


J 


If  we  apply  7 „ to  the  vector  R.  we  get  a column  matrix  of 

^ IJ  13. 

4-  a 

vectors : 


= V ^ R. 
la  fB  la 

**  a 


where  we  have  made  use  of  Equations  (3-3)and(3-5),  If  we  apply 


^7  g to  a row  matrix  of  vectors  we  get 


••• 

pB  la  *^pB  lal  laZ  lai 

^ a ^ a 


V V 

‘Vb  ial  ia2 

^ a ^ a 


pB  lan 
^ a 


"^B 

Note  that  each  of  the  elements  7 ^ ^iak  actually  a column  matrix 

Rt 

of  vectors,  and  hence  ^ is  an  n * n matrix  of  vectors.  The 

la 


transpose  of  this  m^rix  is 


.B  la 
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F |B  ' 
fB  ial 
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We  now  return  to  Equation  (96)  and  take  its  transpose: 
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;B  ia 

^ak  a ,B 
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(4-101) 


Taking  the  transpose  of  this  and  substituting  into  Equation  (95) 
yields 

,t 


' ^B  la  J ^ a 
a 


(4-102) 


Evaluation  of  D 


From  Equations  (90)  and  (102)  we  get  immediately 
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ca 


1 la  a a 


MR  - MR  • uo 
ca  ca 


B 


(4-103) 


E 


m.  $. 

1 la 


7 R 1 
-B  a ' ^ 

’ a 


(4-104) 


where  we  have  made  use  of  Equation(3-17)in  the  form 


a®  = m.$^ 


a . 1 la 

1 


(4-105) 


E 


B 


-B  a 


m.  V n 
1 ^B  xa 


(4-106) 


When  substituting  Equation  (90)  for  R.  into  Equation  (69) 

IcX 
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for  D_2i  we  note  that  in  the  third  term  gets  "caught"  between 

~ ~t  "*R 

R.  on  the  left  and  R.  on  the  right.  In  order  to  pull  i*-'  out  of  the 
la  la  ® , 

summation  over  the  particles,  we  make  use  of  the  following  : 


R.  • UJ  • R.  = • R.  R.  = R.  R.  • 

la  la  la  la  la  la 


Note  that  R.  R.  and  R.  R.  are  triadics.  We  run  into  a similiar 
la  la  xa  la  . 


situation  when  we  substitute  R.  into  Equation  (71)  for  D_,.  In 
~R  gt‘”3 

this  c ase,  <x)  gets  caught  between  R.  on  the  left  and  on  the 

la  la 


right.  In  this  case,  we  use 


~B  . , b'  ^ j-  b'  . JB 
la  la 


When  substituting  Equation  (102)  for  into  Equations  (70) 


and  (72)  to  get  ^33’  respectively,  we  make  use  of 


ia  ia 


The  final  result  is  that  D is  given  as  follows 


”..vm  VR  ..'M 


I /l-  , '''\'  1 '■  ■' 


^ . kR’''.')’'-"  • ^ ,',’1  ■ 


* See  Appendix  C 
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Evaluation  of  X 


We  can  now  get  explicit  expressions  for  X by  using  this  D in 
Equation  (78).  We  first  notice  that  the  first  element  of  X,  name- 
ly X"*  is  zero.  The  second  element  of  X,  namely  , was 

'"a 

simplified  as  shown  in  Equation  (84)  because  of  Equation  (81). 

We  can  explicitly  verify  Equation  (81)  by  use  of  D of  Equation  (110). 
When  we  do  this,  we  must  make  use  of  the  following 


. y . SB  ^ 


m.R.  R.  • a' 
1 la  la 


= w • 


m.R.  R. 

1 la  la 


s®*  E m.T®  • *b‘  5®  = 


Equation  (111)  follows  from  writing  I in  the  form 

4^  \ ^ r -4  4^ 

I = m.  (R.  • R.  ) E - R.  R. 

a . 1 la  la  la  la 

1 


and  Equation  (IIZ)  follows  because  is  an  n x n skew- 

^ la  la 

symmetric  matrix  of  vectors.  It  then  follows  that 


X = -MR  -V  = V X P 
ca  a a 

(1) 


=1*  See  Appendix  C 


J 


Rather  than  writing  out  all  of  X . „ , it  is  convenient  to  just 

i 

th  — ^ th 

write  out  the  j element,  X „ . This  j element  is  obtained  by 

^a- 

th  • J3 

just  using  the  j element  in  the  column  naatrices  a and  and 

th  — a.  jc 

using  just  the  j row  of  the  matrix  023’  Thus,  from  Equation 

(78)  or  (84)  we  have 


"♦R  — -♦R  — -R 

-^!j  • -a  - (°3z)j  • " - <°33)j.  «a 


where 


= E nai.®. 

aj  ^ 1 laj 


(D^^).  _ E • ■ r'^ 

J - i 1 laj  la 


(D  ).  = E ni.t®.  • 

33'j.  ^ 1 laj  la 


When  we  take  the  appropriate  rows  out  of  D as  given  in  Equation 
(110)  we  find 


t ^ I t J a a aj  a 


n.,R„  • • 


-f  i ia  iaji 


1 laj  la 


Em.  -3''  •Em.T.'*. 

. 1 ^^B  lajy  la  ^a  1 laj  la  ^a 


Now  notice  that  the  terni  which  is  quadratic  in  can  be  written 


-S'  ■ • R‘  • S®  = -t'*’  -E'"  1^”  • r' 

1 laj  la  2 i laj  i 


t + R.  -T.^.  • 
la  la  laj 


(4-120) 


and  similiarly  for  the  term  which  is  quadratic  in  ^ •’ 

-f'E-./'t' » 'i  • 5®  = - i-  S®‘E».  \h  V 

^a  A-'  la  ^ a 2 ^ a i laj^  la 


ia  (^^B  iajj  «a 


(4-121) 


In  both  Equations  (120)  and  (121),  we  replaced  a term  by  its 
symmetric  part  because  the  skew- symmetric  part  drops  out  of 
the  quadratic  form  expression. 

Next  we  make  use  of  Equations{3 -3)  and  (98)  and  write 


^ laj  VB  I a fB 
a ^ a \ ^ aj 


- 7„R. 

' ^ aj  a 


aj 


(4-122) 


where  we  have  assumed  that  R.  and  its  first  two  partial  dcriva- 

la 

tives  are  continuous  so  that  the  order  of  partial  differentiation 
can  be  interchanged.  Similiarly  we  get 


aj  Y ' \i  ' 

aj 


— (MR  ) (4-123) 

^ B 

aj 


-f.l- 


a 


7 

^ a 


aj 


w 

aj 


a a 


aj 


(4-124) 


Making  use  of  Equations  (120)  to  (124),  X „ now  becomes 

4 ■ 

aj 


where  T is  given  in  Equation(3-26)  or(3-29).Thus,  we  have  now 
established  that 


X 


T 


(4-126) 


as  already  indicated  in  Equation  (88). 
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Evaluation  of  Y 


To  get  Y we  form  the  transpose  of  D of  Equation  (110)  and 


-^t 


dot  multiply  it  into  a . Since  the  first  column  of  D is  zero,  and 


none  of  the  elements  of  D involve  v 


it  follows  that  Y does  not 
=t 


involve  v . Now  each  element  of  the  second  column  of  D has  a 
term  which  ends  in  a)  , and  this  term  drops  out  when  we  multiply 
it  into  a)  ; therefore,  Y does  not  involve  all  of  D just  as  we 
found  earlier  that  X does  not  involve  all  of  D. 

When  we  multiply  into  a , the  result  is 


Y-»  = uj 

V 

a 


B 


MR 


-*B 

U) 


ca 


+ jB'4,  B*' 

a ^ a ^ a I ^ 


:B 


Y = + 24® 

:^B  a ^ a 


IB 


A 


Y ^ • ^m.T®  • R.  ■ uu®  + 2^®  -^2  4 

113  13  1 1 1 ^ 


B ^B  lB 
11a  la  a 


(4-127) 


‘E-.r*  -T®  •S’t  r 

. 1 la  la  ^a  ^ 1 la  I ^B  la  P 

^ ^ \^a  / 


(4-128) 

(4-129) 


y xB  . /;By  xB'^\  :B 

(^a  ^B^aly 

^ a ' 


In  order  to  separate  from  the  particle  dependent  factors  in 
a ^ 

Y , we  write  out  just  the  j term  as  follows: 

^ a 


aj 


-»B 

U)  • 

E 

m . 9.  . ' 

* R . 

. 4 2^ 

i 

1 laj 

la 

E 

m.t^'  . . 

a 

i 

1 laj 

la  j ^ a 

tB  xB  ;B 
1.  ?.  . • 9.  t 
1 laj  la  a 


(4-130) 
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V.  SUMMARY  OF  DYNAMICS  EQUATIONS 
FOR  A DEFORMABLE  BODY 


In  the  previous  section,  we  have  used  the  transformation 
operator  formalism  to  generate  a set  of  momentum  formulation 
equations  plus  a set  of  velocity  formulation  equations.  The 
momentum  formulation  equations  are  the  following: 


— — n 

-n 

— — 

•4 

p 

0 

F 

u 

-»  -* 

-♦ 

H 

+ 

V X P 

= 

L 

a 

a 

a 

. B 

, B 

Sa 

^;b 

k 

a 

— — 

^ __ 

(5-1) 


and 


1 

4* 

ME 

MR*^ 

ca 

a 

V 

a 

H 

a 

= 

MR 

ca 

4* 

I 

a 

a 

• 

B 

Sa 

a 

e® 

a 

a 

1 

ffl  "J 

1 

The  velocity  formulation  equations  are 


ME 

MR 

ca 

a 

1 

Y-+ 

V 

a 
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TU, 

1 

MR 

ca 
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I 

a 

a 
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= 

-» 
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(5-2) 


(5-3) 
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These  dynamics  equations  form  a complete  set  (in  fact,  two 

complete  sets)  because  they  allow  for  the  determination  (say,  by 

numerical  integration)  of  the  6 external  and  n internal  degree  of 

freedom  motion.  Of  course,  we  must  also  add  some  kinematic 

equations,  but  these  equations  are  straightforward. 

The  dynamics  equations  are  completely  general  because  we 

have  not  introduced  any  constitutive  equations  to  characterize  a 

particular  material.  The  method  of  specializing  these  equations 

is  to  specify  for  i = 1 to  N.  If  we  specify  via  the  eigen- 
la  lag 

vectors  of  the  linearized  dynamics  equations,  then  the  are  the 

modal  coordinates,  and  we  get  the  dynamic  equations  of  Bodley 
Z8  38 

and  Park  ’ ; if  we  linearize  our  dynamics  equations  about 

8 7 

zero  velocity,  we  get  the  linearized  vibration  equations.  In 

BI 

such  a case,  we  would  also  express  k (see  Equations(  4-56)  to 
(4-58)) 

in  terms  of  a stiffness  matrix  and  perhaps  in  terms  of  a 
damping  matrix.  In  this  case,  i.  could  be  considered  as  con- 

1 3^ 

stant  with  respect  to  frame  B (which  has  inertial  angular  velocity 

”*B  R 

^ ),  and  hence,  all  the  partial  derivatives  of  would  be  zero. 

la 

Our  equations  are  really  much  more  useful  than  to  just  use 
them  in  a linearized  form.  As  a matter  of  fact,  the  equations 
can  easily  be  used  in  their  full  generality  in  the  case  of  a col- 
lection of  rigid  bodies.  In  this  case,  the  internal  coordinates 

a 

can  be  used  to  represent  the  relative  coordinates  (rotational  and/ 
or  translational)  between  the  rigid  bodies;  frame  B can  be  fixed 
in  one  of  the  bodies,  say  in  the  main  body;  k now  represents 

3 

interbody  torques  and/or  forces  due  to  actuators  (interbody 

torques  and  forces  which  are  due  to  constraints  drop  out).  The 

resulting  equations  are  given  in  Reference  (69). 

We  can,  of  course,  also  apply  our  equations  to  a collection 

B BI 

of  flexible  bodies  by  introducing  appropriate  and  k . This 


- ( ) 6 - 


subject  will  be  addressed  in  a future  report. 

All  the  dynamics  equations  of  References  (1)  to  (54)  have  the 
general  form  of  Equations(  5-1)  and(  5-Z),  or  of  Equation  (5-3). 
Evidently,  all  of  these  equations  are  straightforward  consequen- 
ces of  Newton's  law  for  a particle:  all  we  need  to  do  is  to  lin- 
early transform  from  particle  velocities,  momenta,  and  forces 
to  system  velocities,  momenta,  and  forces,  respectively.  In 
order  to  obtain  the  extra  term  X or  Y,  we  need  to  evaluate  the 
time  derivative  of  the  transformation  operator  A;  but  this  too 
is  an  entirely  algebraic  process. 

The  equivalence  of  various  forms  of  velocity  equations  has 

88 

been  pointed  out  by  Likins  . In  the  present  work  we  have  ex- 
tended this  equivalence  to  the  momentum  equations  and  velocity 
equations  for  an  arbitrary  deformable  body. 
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VI.  CONCLUDING  COMMENTS 


We  have  seen  that  it  is  possible  to  obtain  the  exact  dynamics 
equations  for  a deformable  body  by  using  purely  vectorial  tech- 
niques. The  result  is  the  same  as  that  obtainable  from  Lagrange's 
or  Hamilton's  equations,  but  we  do  not  have  to  perform  partial 
derivatives  of  the  kinetic  energy  (except,  as  we  did,  to  show  that 
the  results  are  equivalent).  Constraint  forces  (and  torques)  drop 
out  in  our  formulation  for  precisely  the  same  reason  that  they  drop 
out  of  Lagrange's  equation.  In  fact,  we  used  essentially  the  same 
procedure  to  derive  our  equations  as  Lagrange  used  to  derive 
Lagrange's  equations. 

Our  approach  is  very  similiar  to  that  of  Kron  and  that  used  in 

"matrix  structural  analysis",  and  in  the  "finite  element"  method. 

However,  there  are  also  some  differences  between  our  method 

and  that  of  others.  A major  difference  is  that  we  use  classical 

vectorial  mechanics  to  get  all  the  results  we  need.  We  introduce 

angular  velocity  vectors  (rather  than  limiting  ourselves  to  time 

derivatives  of  Euler  angles);  in  this  respect,  our  equations  are 

similiar  to  those  obtained  by  Boltzmann  and  Hamel,  sometimes 

H 0 

called  Lagrange  equations  in  quasi-coordinates.  Such  nonholo- 
nomic  velocities  are  also  used  extensively  by  Kane. 

Finally,  we  should  point  oul  that  we  have  developed  both  n 
set  of  momentum  equations  and  a set  of  velocity  equations.  One 
set  may  be  more  useful  than  the  other  set,  depending  on  the 
particular  problem  at  hand.  Unfortunately,  momentum  formu- 
lations are  unfamiliar  to  many  and  are  not  used  much.  Russell, 
and  Vance  and  Sitchin  arc  conspicuous  exceptions. 
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left  inverse  of  transformation  operator  A 
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■ th  ‘ 

i element 
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particle  i 

total  angular  momentum,  about  point  a, 
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, BE 
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external  part  of  k 

a 

45  (4  - 56) 
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internal  part  of  k^ 

45(4-56) 
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primitive  system  force;  column  matrix 
of  N vectors,  with  f^  as  the  i'''^  element 

33(4-1)  and 
41(4-41) 

K 
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particle  i 
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total  moment  or  torque,  about  point  a, 
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mass,  on  system  of  N particles 

19(2-50) 
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total  external  moment  or  torque,  about 
point  a,  on  system  of  N particles 

19(2-46) 

m . 
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mass  of  particle  i 

9(2-2) 

M 

total  mass  of  system  of  N particles 

9(2-3) 
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n 

number  of  internal  degrees  of  freedom 
of  system  of  N particles 

25 

N 

number  of  particles  in  deformable  body 

9 

o 

zero  dyadic 

48(4-67) 

linear  momentum  of  particle  i 

10(2-5) 

~p 

total  linear  momentum  of  system  of 
N particles 

10(2-6) 

r 

a 

position  vector  to  point  a from  inertial 
reference  origin 

14(2-23) 

r 
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position  vector  to  center  of  mass 

9(2-2) 

r. 
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position  vector  to  particle  i from  inertial 
reference  origin 

9(2-1) 

r. 

1 

inertial  time  derivative  of  r^ 

9(2-1) 
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ca 

position  vector  to  center  of  mass  from 
point  a 

14(2-24) 
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ca 

skew  - symmetric  dyadic  formed  from 
vector  R 

ca 

22(2-63) 
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ca 

dyadic  transpose  of  R 

c a 

22(2-t>2) 

T. 

la 

position  vector  to  particle  i from  point  a 

14 (2-23) 

la 

inertial  time  derivative  of  R. 

la 

14(2-26) 
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R.  ■ R. 
la  la 


16(2-33) 


skew  --  symmetric  dyadic  formed  from 
vector  R. 


dyadic  transpose  of  R.. 


21(2-60) 


time  derivative  of  R.  with  respect  to 

la 

frame  B which  has  inertial  angular 
\'elocity 


21(2-58) 


position  vector  to  particle  i from  center 
of  mass 


14  (2-27) 


kinetic  energy  of  particle  i 


10(2-8) 


total  kinetic  energy  of  system  of  N particles  10(2-9) 


kinetic  energy  function  in  terms  of  v 


kinetic  energy  function  in  terms  of  v 


kinetic  energy  function  in  terms  of  v 
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alternate  notation  for  R. 
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la  la 
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velocity  of  point  a 

14(2-26) 

-►2 

16(2-33) 
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a 

V • V 

a a 
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velocity  of  center  of  mass 

10(2-4) 

V. 
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velocity  of  particle  i 

9(2-1) 
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inertial  time  derivative  of  vector  V 

21(2-57) 
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— ► 
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time  derivative  of  vector  V with  respect 
to  frame  B which  has  inertial  angular 
velocity  ^ 

21  (2-57) 

X 

a primitive  term  which  is  cjuadratic  in 
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^3(4-1) 

X 
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36  (4  - IH ) and 
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X-. 

V 

first  part  of  the  three-p;irt  column- 

51(4-84) 
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decomposition  of  X 

second  part  of  the  three-part  column- 
decomposition  of  X 

51(4-84) 

third  part  of  the  three-part  column- 

51(4-84) 
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decomposition  of  X 
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aj 

element  of  Xig 
'’a 

51  (4-115) 
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a primitive  term  which  is  quadratic  in 

34  (4-3) 

CJ  or  in  G 
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a transformed  term  which  is  quadratic 
in  ■'  and  G 

36  (4-21)  and 
50(4-79) 

Y . 

first  part  of  the  three-part  column- 

52(4-86) 
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decomposition  of  Y 

YoJB 

second  part  of  the  three-part  column- 
decomposition  of  Y 

52  (4-86) 

third  part  of  the  three-part  column- 
decomposition  of  Y 

52(4-86) 

element  of  ^ 
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coefficient  of  ^ in  expansion  of  internal 
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generalized  momentum  g in  terms  of 
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Va>  , and 
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transpose  of  B 
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special  case  of  B when  a is  the  center 
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3^ 

^ B 

coefficient  of  in  expansion  of 

29(3-19)  and 
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internal  generalized  momentum  g 

— ► ~^B  * B ^ 

in  terms  of  v , , and  ^ 

a a 
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31(3-33) 

a 

R 

transpose  of  0 

a 

29(3-20) 

3^^ 

C 
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special  case  of  3 , when  a is  the  center 
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of  mass  c 

47(4-62) 

a 

coefficient  of  4 in  expansion  of  internal 

^ B 

generalized  mon'ientum  g in  terms  of 

-►  -♦B  , f B ^ 

V , a'  , and  t, 

a a 
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■Q 

special  case  of  y , when  a is  the  center 
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of  mass  c 

47(4-62) 

primitive  system  mass;  diagonal  NXN 
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transformed  system  mass 
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inverse  of  P 
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V 

inverse  of  p 
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th 

j internal  generalized  coordinate  relative 
to  point  a and  frame  B 
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column  matrix  of  n scalars,  with  £^.  as 
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j elcnient 
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transpose  of  ^ (i.  e.  , row  matrix 

^ th 

of  n scalars,  with  4 . as  j element) 

aj 

special  case  of  , when  a is  the  center 
of  mass  c 


primitive  system  velocity;  column  matrix 

^ th 

of  N vectors,  with  v.  as  the  i element 


transformed  system  velocity;  column 


matrix  whose  elements  are  v , 

, £ B ^ 

and  s 
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. coefficients  of  expansion  of  u.  in 
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terms  of  (equal  to  partial  derivative 

of  R.  with  respect  to  ) 
la  ^ ^aj 
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laj 
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laj 


skew- symmetric  dyadic  formed  from 
vector  $ ! 


B 

iaj 


dyadic  transpose  of  $ . 


B 

iaj 


$ : 


B 


la 


B^ 

ia 


B 


column  matrix  of  n vectors,  with 

. as  element 
laj 

transpose  of  ( i,  e.  , row  matrix 

of  n vectors,  with  '1.  .asi  element) 

laj 

special  case  of  ^ . when  a is  the 
center  of  mass  c 


(equation) 

occurrence 

27(3-8) 

47(4-62) 

33(4-2)  and 
39(4-38) 

35(4-12)  and 
42(4-47) 

25(3-1) 

54(4-91) 

54  (4-93) 
26(3-5) 

26(3-6) 

47(4-63) 
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NOMENCLATURE  - Continued 


j 

I 


Symbol 

Definition 

Page  (equation) 
of  first  occurrence 

la 

column  matrix  of  n dyadics  formed 
from  by  replacing  all  vector 

53(4-91) 

elements  by  the  skew- symmetric 
dyadic  of  these  vectors 

ia 

r 7B 

transpose  ot  ? . 

^ la 

54(4-93) 

$ ® 
la 

time  derivative  of  ^ with  respect  to 

frame  B 

55(4-95) 

angular  velocity  of  frame  B with 
respect  to  inertial  space 

Zl(2-57) 

skew- symmetric  dyadic  formed  from 
vector  ^ 

53(4-89) 

VfB 
^ a 

column  matrix  of  n scalar  partial 
differentiation  operators,  with 

55(4-97) 

S .th  , 

^ as  j element 

SC® 

aj 

0 

zero  vector 

12(2-18) 

APPENDIX 


A.  VECTORS,  DYADICS  AND  TRIADICS 

Our  physical  world  includes  objeits  which  we  like*  to  count 
or  measure.  For  this  purpose  the'  inti’gers  and  ri.’al  numbers  were 
invented.  Originally  these  numbers  did  not  exist  (only  the  counting 
and  measuring  existed);  but  ev<.mtiuilly  these  numbers  were  ac- 
cepted as  existing  mathematically  because  they  t ould  be  used  tor 
counting  and  measuring.  These  numbers  were  used  extensively 
in  the  geometry  arid  physics  ot  our  i-dimensional  world.  Many 
geometric  and  physical  rjuantities  actually  have  .3  real  numbers 
associated  with  them,  and  therefore  these  geometric  and  physical 
quantities  were  eventually  represented  by  3-dimensionaI  vectors. 
Originally,  these  vectors  did  not  exist  (only  the  3 real  numbers 
existed);  but  eventually  these  vectors  were  accepted  as  existing 
ma thetJta ti ca  1 1 y because  they  i ould  be  used  lor  producing  the  3 
real  numbers  associated  with  the  geometrical  and  physical  quanti- 
ties. Siniiliarly,  dyadics  eventually  were  accepted  as  existing 
because  they  could  be  used  to  generate  vectors;  ;ind  triadic  s weri' 
aciepted  as  existing  because  they  could  be  used  to  gencM'ate  dya- 
dics; etc. 

Vectors 

If  X,  y,  and  z are  used  to  symbolii  ally  represent  5 orthogon- 
al  directions  in  our  3-dimensional  space,  then  a vector  V can  be 
written  as 

V = V X t V y I V /.  (A-1) 

X y z 

Note  that  we  ,i  ri'  not  really  adding  V^,  V^_,  and  V^  ; nor  are  we 
adding  x,  y,  and  z.  Efjuation  (A-l)  is  merely  a symbolic  repre- 


- A\- 


sentation  of  the  fact  that  the  geometric  or  physical  quantity  h.  s 

the  3 real  numbers  V , V , and  V in  the  3 directions  x,  y,  ^nd 

X y z 

z,  respectively.  The  operation  of  the  dot  product  is  then  deiiued 
so  that  we  can  produce  real  numbers  from  vectors.  Thus,  for 
the  3 orthogonal  directions  x,  y,  z we  define 


A 

X 

• X = 1 

A . 

X 

'<> 

II 

o 

A 

X 

A 

Z 

A 

A ^ 

A 

II 

A . 

A 

y 

• X = 0 

y 

y 

Z 

N> 

x> 

II 

o 

A , 

Z 

II 

o 

A . 

Z 

A 

Z 

(A-2) 


Then  we  get  the  vector  components  V , V , V as  follows 
° X y z 


V 

- A • 

V = V • 

A 

X 

X 

V 

= ^ • 

V = V • 

.A 

y 

(A-3) 

y 

V 

Z 

A . 

= Z 

V = V • 

Vector  components  can  be  conveniently  put  into  matrix  form. 
Thus,  we  define  and  as  follows 


V 


Q 


V 

y 


V 


V 


t 

Q 


V 

y 


V 


(A-4) 


We  can  similiarly  put  the  3 directions  x,  y,  z in  a matrix  form. 
Thus,  we  define  Q and  Cl  as  follows 


„t  rA  A A 

Q = Lx  y z J 


Using  the  ordinary  rules  for  matrix  multiplication,  Equation  (A-1) 


can  now  be  written  as 


V = Q = qVq 


(A-6) 


Similiarly,  Equation  (A-2)  can  now  be  written  as 


Q • Q = h 


(A-7) 


where  1^  is  the  3x3  unit  matrix.  Equation  (A- 3)  can  now  be 
written  as 


Vq  = Q • V = V • Q 


(A-H) 


One  might  argue  that  matrices  were  not  intended  to  be  used 
with  vectors  as  elements  as  in  Equation  (A-S).  The  only  response 
to  such  an  objection  is  that  one  then  simply  redefines  the  notion 
of  a matrix,  but  uses  the  same  old  matrix  symbols,  just  as  in 
Equation  (A-l)  we  tacitly  redefined  the  notion  of  addition  so  that 
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it  would  also  apply  to  vectors. 
Dyadic  s 


Physical  quantities  which  have  9 real  numbers  associated 
with  them  are  conveniently  represented  by  a dyadic.  Thus,  the 
dyadic  D can  be  written  as 


D = D XX  + D xy  + D 

XX  xy  X2 


+ 

D 

A A 

yx 

+ 

D 

A A 

yy 

+ 

D 

A A 

yz 

yx 

yy 

yz 

+ 

D 

A A 
ZX 

+ 

D 

+ 

D 

A A 
ZZ 

zx 

zy 

zz 

(A-9) 


AAA 

We  denote  the  x,  y,  z components  of  the  dyadic  D by  the  matrix 
Dq  given  by 


D 

D 

D 

XX 

xy 

xz 

D 

D 

D 

yx 

yy 

yz 

D 

D 

D 

zx 

zy 

zz 

(A-10) 


D 


can  now  be  written  as 


^ - q‘  DqQ  (A-11) 

Since  Q * Q*'  = 1^  we  can  solve  this  for  as  follows; 

Dq  = Q • *D  • (A-12) 
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Thus,  whereas  it  takes  one  dot  multiplication  with  Q to  produce 
a set  of  real  nunibers  from  a vector  (see  Equation  (A-B)  ),  it  takes 
two  dot  multiplications  with  Q to  produce  a set  of  real  numbers 
from  a dyadic. 

A dyadic  times  a vector  is  another  vector.  Thus, 


D 


= (Q^DqQ)  • (Q^Aq) 


Q^B 


Q 


(A-13) 


J 


Thus  we  can  write 

D • A = B implies  D^A^  = Bq 


(A-14) 


Note  that  the  side  from  which  D is  multiplied  is  important 
because 


A • ^ = (A^  Q)  • (Q^DqQ) 

= (A^  Dq)  Q 
= q'(d^  Aq) 


(A-15) 


Thus  we  can  write 


3 
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A • D = C 


(A-16) 


( 


I 

1 

\ 


1 


[ 

I 


i 

I 


implies  Aq  = 


The  transpose  or  conjugate  of  D is  the  dyadic  D such  that 


D • A = A • D 


(A-17) 


We  evidently  have 


t t t 

d'  = Q Dq  Q 


(A-18) 


Thus,  the  matrix  of  components  of  D is  the  transpose  of  the 
matrix  of  components  of  D.  If  D equals  D then  D is  called  a 
symmetric  dyadic;  if  D equals  -D  then  D is  called  an  anti-sym- 
metric  or  skew- symmetric  dyadic.  Thus,  D has  the  same 
symmetry  as  its  matrix  of  components. 

We  can  form  a dyadic  by  Juxtapositioning  two  vectors 


C = A B = (Q^Aq)  (Bq  Q) 
= q'(AqB^)Q 


(A-19) 


Thus 


C = AB  implies  Cq  = A^  Bq  (A-20) 


Note  that  in  general,  AB  is  not  the  same  dyadic  as  BA.  Even 
though  we  can  form  a dyadic  out  of  two  vectors,  riot  every 
dyadic  can  be  formed  this  simply.  In  the  general  case,  it  takes 
three  vectors  to  specify  a dyadic.  For  exaniple,  the  dyadic  I) 


I 


* f 
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r 


of  Equation  (A-9  ) can  be  written  as 


^ A A A 

D=xV  iyV  +zV 
X y z 


(A-21) 


whe  re 


V = 

D 

A 

X 

+ 

D 

A 

y 

+ 

D 

X 

XX 

xy 

— > 

V = 

D 

A 

X 

+ 

D 

A 

y 

f 

D, 

y 

yx 

yy 

II 

t> 

D 

A 

X 

+ 

D 

A 

y 

+ 

D 

Z 

zx 

zy 

yz 


(A-22) 


Thus  the  3 vectors  V V , and  V (together  with  the  basis  vec- 

X y z 

tors  X,  and  z)  represent  the  dyadic  D.  Alternatively,  we  can 

represent  D by  W , W , and  W as  follows: 

^ X y z 


■4— > — > . A A 

D=W^-tWy  + Wz 
X y z 


(A-23) 


wh  ere 


w = 

D 

X + D 

A 

y 

+ 

D 

X 

XX 

yx 

zx 

-> 

w = 

D 

.<$  + D 

A 

y 

+ 

D 

y 

xy 

yy 

zy 

w = 

D 

X 1 D 

A 

y 

+ 

D 

z 

xz 

yz 

zz 

(A-2-4) 


For  any  vector  A there  is  a dyadic  A such  that  the  vector 
cross  product  becomes  a dyadic  dot  product.  Thus  lor  any  A 
and  B: 
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It 

ri 


ij 


A X B = A • B (A-25) 

Since  A x B = -B  x A we  also  have 

A X B = -B  - A (A-26) 


If  A , A , A are  the  x,  y,  z components  of  A,  then  A can  be 
X y 2 

.AAA  , / 1 1 

given  in  x,  y,  z components  as  loilows 


A 


A AA,  A AA,  A AA 

==  -A  xy  + A xz  + A yx 

z y 7. 

A AA  A AA  ,a  AA 

-A  yz  - A zx  1 A zy 

X ' y X ^ 


(A-27) 


Thus,  we  can  write 
A = Q^AqQ 


where 


Since  Aq  is  skew- symmetric , so  is  A: 


A = -A 


(A-28) 


(A-29) 


(A-iO) 


Using  Equations  (A-17)  and  (A- 30),  Equations  (A-25)  and  (A-26) 
can  now  be  written  as  follows 


1 
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2<f2 

A[>AO42950 


(A-31) 


-4 

A X B = A • 

-4 

B 

-4 

= B - 

A 

t 

•4 

-4 

= 

• A 

= A • 

B 

The  identity  dyadic  E is  defined  such  that  for  any  vector  A : 

-♦  -»  ■*-*  •* 

A • E = A = E 'A  (A-32) 

E can  be  expressed  in  terms  of  A,  y,  z as  follows 

E = Q^Q  = XX  + ^y  + zz  (A- 33) 

t 

If  we  write  E = Q E^  Q then  E„  = 1,,  the  3x3  unit  matrix.  The 
zero  dyadic  O is  defined  such  that  for  any  vector  A: 

A • O = 0 = O • A (A-34) 

•*-»  f 

If  we  write  O = Q Oq  Q then  Oq  is  the  3x3  zero  matrix. 

We  can  write  the  well-known  vector  triple  product  identity 
in  the  following  form: 

A • B • C = A X (B  X C ) = (A  • C )B  - (A  • B ) C (A-35) 

= [BA  -(B  • A)E]  - c 

From  this  it  follows  that 

A-B=BA-(B-  A)E  (A-36) 

Next  we  write 
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(AxB)  • C=(AxB)kC 

= -C  X (A  X B ) 

= C * (B  X A ) 

-»  -»  -♦  -»  -♦  -* 

= (C  • A )B  - (C  • B )A 

= c • Cab  - b a ] 

= [B  A - A B ] • C 


(A-37) 


In  the  last  line  we  made  use  of  the  fact  that  the  transpose  of 

-*  -♦  -•  -* 

A B is  B A : 


D = A B 


implies  D = B A 


(A-38) 


Hence  we  have 


AxB  =BA-AB 
= A • B - B 


(A-39) 


T riadics 


A triadic  T has  the  general  form 


AAA 

xxz 


+ T ^zx  + T ^zy  + T ^zz 
xzx  xzy  ^ xzz 


T 

T 

^xy  + T 

XXX 

xxy 

xxz 

+ T 

T 

xy^  + T 

xyx 

xyy 

' xyz 

rr'  A A ^ I 

+ T yxx  + , . , 
yxx 


+ . . . 


I AAA  , 

+ T zxx  + . . . 
zxx 


(A-40) 


A special  triadi 


c is"*^ 


defined  as  follows 


AAA  , AAA  I AAA 

-zxy  + yxz  + zyx 

AA/»  AAA  , AAA 
-xyz  - yzx  + xzy 


(A-41) 


By  direct  expansion  we  find 

- X ~ 

A = A 


(A-42) 


Z ' B = B 


(A-43) 


A • Z • B = A X B 


(A-44) 


Thus  the  triadic  Z permits  us  to  give  equal  treatment  to  A and 
B in  A X B ; i.  e.  , we  can  let  both  of  them  be  vectors,  rather 
than  converting  one  of  them  to  a dyadic. 

Let  be  the  Levi-Civita  symbols  defined  as  follows 

(+1,  if  (i,  j,  k)  is  an  even 

permutation  of  (1,  2,  3) 


Also  let 


-1,  if  (i,  j,  k)  is  an  odd 

permutation  of  (1,  2,  3) 

0,  otherwise  (if  two  or  more  of 
i,  j,  or  k are  equal) 


-All- 


(A-45) 


(A-46) 


Z can  now  be  expressed  as 


i'- 


J .J  V 

- E L i: 

i = 1 j = 1 k = 1 


, AAA 

f . e. e .e, 
ijk  1 j k 


(A-47) 


Note  that  E is  similiarly  related  to  the  Kronecker  symbol  6. 


ij 


where 


^iJ  = 


1,  if  i = j 


0,  otherwise  (if  i ^ j) 


(A-48) 


Evidently,  E can  be  written  as 


E=  E E 


i = 1 j = 1 


tj  t J 


(A-49) 


Time  Derivatives 

We  will  now  obtain  expressions  for  the  inertial  time  deri- 
vatives of  vectors  and  dyadics  in  terms  of  the  time  derivatives 
with  respect  to  a rotating  frame.  Consider  a frame  denoted  as 
"frame  B"  with  rectangular,  orthogonal,  right-handed  unit  vec- 
tors x„,  y,.,,  z„.  Then  the  vector  ^ can  be  written  as 
B D rS 


V = Q V 


(A-50) 


B 


where 
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B 


Q 


B 


B 


B 


(A-51) 


and 


B 


‘B 


(A-52) 


'B 


If  the  inertial  angular  velocity  of  frame  B is  U)  , then 


A **B  „ A ~ B A A ~B 

’^B  " Xb  = “)  • Xg  = Xg  • u) 


(A-53) 


and  similiarly  for  y and  z„.  Hence 

ti  B 


=pB  ^b]  >*B  ^b] 


= • Qb 


Taking  the  inertial  time  derivative  of  V now  yields 


7 

~B  t t • 

= ou  • Q V + O V 

bX  ® Qb 

= uo®  • V + V 


(A-55) 
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where  V is  the  time  derivative  of  V with  respect  to  frame  B (i.  e. 
the  time  derivative  obtained  if  frame  B is  considered  fixed). 

Thus,  we  have 

^ S I A r L \ 

V = V + uu  X V (A-56) 


For  a dyadic  D 


we  get  similiarly 


° ^B  *^B  ^B  ^ 


UD  • D + D • uu-“  + D 


(A-58) 


where  we  have  used  the  transpose  of  Equation  (A-54): 

>t 


^B  - ^B 


u: 


rB 


(A-59) 


and  we  have  defined  D as  the  time  derivative  of  D with  respect 
to  frame  B.  Evidently,  the  time  derivative  of  the  identity 
dyadic  is  zero 


(A-60) 


•♦B 

Note  that  Equation  (A- 55)  has  one  term  in  u)  , and 

O 

Equation  (A-58)  has  two  terms  in  uu  Clearly,  for  the  inertial 

•♦B 

time  derivative  of  a triadic  we  will  get  three  terms  in  (d  . 


-A14- 


B.  MATRICES  OF  DYADICS,  VECTORS,  AND  SCALARS 


It  is  often  convenient  to  mix  matrices  of  scalars,  vectors, 
and  dyadic s.  For  example.  Equation  (A -22)  can  be  written  in 
matrix  form  as 


V - DqQ 


(B-1) 


where  Q and  Dq  are  as  in  Equation  (A-5)  and  (A-10),  and  V is 
given  by 


V = 


(B-2) 


Now  let  be  defined  by 


V 


t 


(B-3) 


Then  direct  expansion  shows  that 

V^=Q^Dq  (B-4) 

But  this  result  can  also  be  obtained  directly  from  Equation  (B-1) 
by  just  formally  taking  transposes.  The  validity  of  such  formal 
manipulations  must  be  justified,  but  the  justification  merely 
involves  direct  expansion. 

Suppose  V and  W are  column  matrices  of  two  vectors  and 
two  scalars 
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^1 

-> 

^2 

; W = 

-* 

W 

2 

(B-5) 

W3 

74 

.^4 

1 

Here  V^,  V21  vectors  and  V^,  V^,  W^,  are 

scalars.  Next  suppose  we  have  four  equations  as  follows 


= ^11  • 

+ A^2 

W2  + Ai3W3  + .?14W4 

-# 

-4  -»  -4 

V = A 

21 

^1  ^ "^22 

• ^2  + ''23^3  4 ^24  "^4 

•4 

^3"  '^31  ‘ 

Wi  + A32 

+ A33  W3  + A34  W4 

V.  = A,,  • 
4 41 

Wi  4 A42 

+ A43W3  + A44W4 

(B-6) 


It  is  evident  that  we  can  put  these  equations  in  matrix  form  as 
follows: 


* ^ 

44 

^11 

44 

'^12 

Ai3 

4 

^14 

-> 

^2 

= 

44 

^21 

44 

^22 

■^23 

■^24 

• 

^2 

^3 

A31 

'^32 

A33 

A34 

^3 

^4 

4 

^41 

^^42 

^43 

A. . 
44 

74. 

4-4 

Since  A 

II 

• A* 

All  ’ 

44 

^12 

-» 

• W 

2 

-+ 

= W 
2 

44 

• A 

(B-7) 


t 

12 


-\i2- 


r 


also  write  this  as 


\ ^3  ^4j 


W3  W^J 


^11 

^21 

A31 

■^32 

"^42 

5.3 

^23 

A33 

A43 

A.4 

^^24 

A34 

^44 

(B-8) 


Thus  Equations  (B-7)  and  (B-8)  can  be  written  as 


V = A • W 


(B-9) 


(B-10) 

where  A and  A^  are  the  matrices  in  Equations  (B-7)  and  (B-8), 
respectively.  Note  that  the  dyadic  elements  of  A were  changed 
to  their  transposes  when  A was  changed  to  A*'. 

In  order  to  have  a uniform  procedure  for  forming  transposes 
of  matrices  whose  elements  may  be  dyadics,  vectors,  or  scalars, 
we  adopt  the  following  convention:  the  transpose  of  a scalar  and 

of  a vector  are  equal  to  the  scalar,  and  vector,  respectively: 


V = V 


(scalar  s) 
(vector  V ) 


(B-11) 
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Consequently,  is  formed  from  A by  taking  the  element  in  the 
i^^  row  and  column  of  A and  putting  the  transpose  of  this 
element  in  the  row  and  i^^  column  of  A^.  Note  that  the  same 
procedure  applies  to  forming  V^'  and  from  V and  W,  respec- 
tively. 

Suppose  X is  a column  matrix  of  n vectors 


(B-12) 


We  now  define  X by 


(B-13) 


Note  that  X^  was  defined  in  Appendix  A for  any  vector  X^.  Now 
we  have  defined  X for  any  column  matrix  of  vectors.  Next,  we 
note  that  our  above  convention  on  transposes  requires  that  the 
transpose  of  X must  be 
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(B-14) 


Thus,  is  the  negative  of  what  one  obtains  by  first  taking  the 
row  matrix  of  vectors  and  then  replacing  each  of  the  row 
elements  X^  by  X^.  The  potential  confusion  here  due  to  this 
change  in  sign  can  be  eliminating  by  simply  remembering  that 
we  start  with  a column  matrix  of  vectors,  then  perform  the  "tilde 
operation",  and  only  then  do  we  take  the  transpose.  In  other 
words,  the  "tilde  operation"  is  only  defined  for  a single  vector 
or  for  a column  matrix  of  vectors  (but  not  for  a row  matrix  of 
vectors). 

-4  -♦  ~ "♦  t 

From  Equation  (A- 31)  we  have  AxB=A’B=B  *A.  A 
generalization  of  this  is  the  following 


AxX  = A • X = A • 


(B-15) 


Also 


-B5- 


A • X^  = A • 

[^1 

^2 

•4 

X 

n 

- 

= A ' 

[^1 

^2 

• * ' 

X 

n 

^2 

X* 

n 

1 

• A 

— f 

= X • A 


A further 


generalization  is 


Y X = Y 


(B-17) 


= -Y  • 


X 


t 


A somewhat  similiar  situation  occurs  if  we  have  a vector 
V given  by 


-B6- 


Am. 


= 2^  X.s.  = s^X  = xS 
1 1 


(B-18) 


where  the  s.  are  scalars  and  s and  s are  column  and  row  matrices  of 
1 ^ 

these  scalars.  We  now  form  V as  follows 


= ^ X.s.  = s^5 
1 1 


"l  "2 


= Pl 


(B-19) 


Equations  (B-16),  (B-17),  and  (B-19)  have  negative  signs  because 
of  the  fact  that  the  "tilde  operation"  and  the  transpose  operation 
do  not  commute  (in  fact,  they  anti-commute). 

Time  derivatives  of  column  matrices  of  vectors  and  row 


-B7- 


matrices  of  vectors  are  easily  obtained  from  Equation  (A-55)  or 
(A-56).  We  find 


and 


B 

X = X 

+ uT® 

• X 

B 

~B^^ 

U) 

= X 

+ X • 

x'  --  x' 

• x' 

J- 

+ x^ 

(B-20) 


(B-2I) 


where  UJ  is  the  inertial  angular  velocity  of  a particular  reference 

B 

frame  (called  frame  B),  and  X is  the  column  matrix  of  vectors 


, th  9 

whose  j element  is  X (if  X is  as  given  in  Equation  (B-12)  ), 

J 

which  is  the  time  derivative  of  X^  with  respect  to  this  particu- 
lar reference  frame. 


1 
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c.  DETAILED  EVALUATION  OF  D AND  p. 

We  will  now  apply  the  formalism  of  the  previous  two  appen- 
dices and  derive  the  results  given  in  Section  IV  D.  Before  doing 
this  we  recall  the  following  definition  of  D given  in  Section  IV  C. 


D = 


o 

O 

0 

MR 

ca 

^22 

°23 

a 

°32 

°33 

(C-1) 


where 


MR 

ca 

= 

E 

m. 

1 

^ia 

i 

E 

m. 

a 

i 

1 

la 

E 

m. 

R.  • 

R.^ 

22 

i 

1 

la 

la 

D,_ 

E 

m. 

. 

R.^ 

32 

i 

1 

la 

la 

Do, 

E 

m. 

R.  • 

23 

i 

1 

la 

la 

D33 

= 

E 

m. 

1 

• 

la 

la 

(C-2) 

(C-3) 

(C-4) 

(C-5) 

(C-6) 

(C-7) 


In  addition  there  are  the  relationships 


•*  — — t •*"*t 

1 = + D,,  = I 

a 22  22  a 


(C-8) 


-CT- 


(C-9) 


P®  = D,,  + D_* 

a 32  23 

■y®  = D + D *■  = y® 
a 33  ^33 


t 


whe  ] 


m.R. 


~t  ♦♦t 
R = I 


a 

1 la 

1 

la 

a 

Sm.fB  ■ 

R.^  = 

m.R. 

a 

1 la 

1 

la 

. 1 la 

1 

la 

y^  = 

y®' 

a 

^ 1 la 

1 

la 

a 

(C-10) 

(C-11) 

(C-12) 

(C-13) 


In  these  equations  the  summation  is  over  all  particles  (from  i = 1 
til 

to  N)  and  the  i particle  is  at  position  R.  relative  to  the  arbi- 

^ la 

trary  reference  point  a.  R.  is  the  inertial  time  derivative  of 

la 

R.  and  can  be  expressed  as  follows 


R.  = R.  + X R, 
la  la  la 


(C-14) 


+ 0)® 
la 


8 

where  R.  is  the  time  derivative  of  R.  with  respect  to  frame  B 
ia  •♦B 

which  has  angular  velocity  u)  with  respect  to  inertial  space. 


^ is  a row  matrix  of  scalars,  each  of  these  scalars  being  the 
time  derivative  of  an  internal  generalized  coordinate 


(C-15) 
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g 

$.  is  a column  matrix  of  vectors 
la 


where 


(C-16) 


a R. 
la 


Thus,  is  a column  matrix  of  the  partial  derivatives  of  R.  with 
la  ^ la 

respect  to  the  internal  generalized  coordinates.  This  fact  can  be 
brought  out  more  clearly  by  introducing  the  operator  {7  as 
follows:  ’a 


uB  = 


Equations  (C-16)  to  (C-18)  can  now  be  combined  as  follows: 


la  VB  la 
^ a 


(C-19) 


We  define  the  operation 
follows 

of  7 on  a 
<! 

row  matrix  of  vectors  ; 

i 

' i 

= 

la 

^ a 

['Si 

... 

idiC, 

] I 

lanj  j 

j 

(C-20)  j 

= 

k.'i. 

^ a 

V 

'1b  ia2 
^ a 

■ • i'b^J  : 

^ a 

Note  that  V -a  ^ ■ is  an  n x n matrix  of  vectors  and  the  element  in 

^ O 1 cL 

the  row  and  k''^  column  is 


(C-22) 


L 


t 1*1  til 

and  the  element  in  the  j row  and  k column  is 


i^R. 
ak  ’ aj 


(C-23) 


Comparing  Equations  (C-21)  and  (C-23)  we  note  that  if  we  f~an 
interchange  the  order  of  the  partial  derivatives  then  the  matrix 
^7  ^ is  symmetric.  The  interchange  of  the  order  of  the  par- 

•S  . . -» 

tial  derivatives  is  permissible  if  and  its  first  two  partial 

derivatives  are  continuous  (i.  e.  are  continuous  functions  of  the 

internal  generalized  coordinates). 

B 

Since  is  a column  matrix  of  vectors  we  have 
la 


= f B + 1®  • 

la  la  la 


(C-24) 
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B 

where  is  the  time  derivative  of  in  frame  B.  But  we  can 
la  la 

write 


la 


k = 1 3 4 


ia  iB 
B ^ak 
ak 


(C-25) 


This  expression  can  be  written  in  terms  of  if  we  first  take 

. ^ a 

the  transpose 

,t 


B 


fB;.  SiB  llia 

k = 1 34 


, V 

ak  i ^B  ^ a .B  ia 


(C-26) 


ak 


Taking  the  transpose  of  this  and  combining  with  Equation  (C-24) 
yields 


la  ( ^B  la  1 ^a  la 


(C-27) 


We  will  use  this  equation  in  Equations  (C-2)  to  (C-13),  as 
required. 

In  order  to  evaluate  Equations  (C-2)  to  (C-10)  we  also  need 
an  expression  for  R.  . Now  from  Equations  (B-18)  and  (B-19) 
we  have 


(C-28) 


where 


-(  6- 


(C-29) 


From  Equation  (A -39)  we  have 


Su®  X R.  = u)®  • R.  - R.  • uj® 

la  la  la 


Combining  Equations  (C-14),  (C-28),  and  (C-30)  now  yields 

R.  = • R.  - R.  • 

la  a la  la  la 


We  are  now  ready  to  evaluate  Equations  (C-2)  to  (C-10). 

Evaluation  of  MR 
ca 

Multiplying  R.  by  m.  and  summing  over  all  the  particles 

1 cL  X 

yields  immediately 


R =2  ni.R.  = » 

ca  ‘—r  1 la  a 


4-  • 


+ 0)“  • MR  - MR  • U) 
a ca  ca 


where  we  have  made  use  of  the  expression 


a®  = m J 

a “-r  1 


(C-33) 
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from  which  it  follows  that 


m.  ? . 

1 la 


(C-34) 


Evaluation  of  D 


22 


From  Equation  (C-4)  and  (C-31); 
D = E "'iRu  • Ria  = « a S >"i*f 


lU  + 

la  la 


•'*1  -]^m.R. 

a 1 1 


~ t 

fii  *R 
la  la 


(C-35) 


where  we  have  introduced  1 from  Equation  (C-H)-  The  last 

3. 

term  on  the  right  hand  side  of  Equation  (C-35)  is  not  in  an  accept- 
able form  since  we  want  to  separate  the  terms  involving  the  in- 
dividual particles  (i.  e.  terms  involving  the  particle  subscripts) 
from  terms  like  ai®  which  do  not  involve  the  individual  particles. 
In  order  to  do  this  we  expand  as  follows,  making  use  of 
Equation  (A-36). 


R.  • uu®  • r!^  = R.  • 

la  la  la 


-*B\  ^ 

-»B 

(R. 

• tio^)  E - 

R. 

\ la 

/ 

la 

(C-36) 


/•*  “»R\  ~ ~ ^ 

= (R.  • UJ  R.  - R.  • R. 

\ la  / la  la  la 


B 


But  the  last  dyadic  on  the  right  is  zero  because  if  we  dot  multiply 
it  into  an  arbitrary  vector  V , we  get  R^^  * R.^  • V/  - 0 


la 


because  R.  x R.  = 0 . Consequently  we  have 
la  la 


R. 

la 


u)“  • • R.  R.  = R.  R. 

la  la  la  la  la 


U) 


B 


(C-37) 
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Note  that  R.  R.  and  R.  R.  are  triadics.  D _ can  now  be 
la  la  la  la 

expressed  as 

D = m R • • rI"  + *"*?  - m.R. 

^2Z  ii3i  ^ia  ^a  4-»  1 la  la  a ^ i i 

1 1 ^ 


R. 
la  la 


(C-38) 


Evaluation  of  D 


23 


From  Equation  (C-6)  and  (C-31): 

23  A-'  i ia  la  ^ a 4-<  i la  la  a “ i i 


w • w. 

la 


(C-39) 


B 

where  we  have  introduced  3 from  Equation  (C-12).  Again,  the 

^ -*B 

last  term  or  *he  right  must  be  manipulated  so  that  ^ can  be 
factored  out.  In  order  to  this  we  use  the  following  from 
Equation  (B-16) 

~B  . jB*^^  l-B^  . ;JB  (C-40) 

ia  ia 


Analogously 


U) 


B 


bb‘  , bb' 
a a 


U) 


►B 


(C-41) 


From  Equations  (C-40)  to  (C-42)  we  now  have 


rB 


R. 

1 la 


~B  . B 
uu  • $. 

la 


t 


«'-E 

a 

1 

= -^mi 


m.R. 

1 la 


la 


R.  • uu 

la 


la 

B 


uo 


(C-43) 


0^2  can  now  be  expressed  as 


^23  4-'  i ia  ia  ^a  4^  i ia  la  ^ i la  la 


(C-44) 


Explicit  Verification  of  Equation (4- 81) 

We  can  use  Equation(4-81)  of  Section  IV  C to  check  the 


correctness  of  the  above  explicit  expressions  for  ^^2  ^23' 

When  we  dot  multiply  of  Equation  (C-38)  into  u)  we  get  three 


terms,  and  the  second  and  third  terms  cancel,  as  can  be  seen 
from  the  following 


(JO®  • 3®  = uj® 

^ i 


= U)' 


i 

= ^ m.  (r.  R 

1 ' la 


;rt 

?B 

• 

R.  • 

uu 

La 

la 

-4 

•4  ^ 

R. 

• R'  1 

1 E - 

la 

la' 

•4 

-» 

-*B 

R. 

R.  • 

U) 

. la 

la 

► 

-*B\ 

~B 

ia 

til  j • 

(JO 

la 


ou 


B 


-»B 

= ou 


rB 


• H m.R 

i 


....  R.  • u) 
1 la  la 


B 


R.  R.  • u) 
la  la 


B 


(C-45) 
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mm 


— “»B 

Hence,  ^ reduces  to 


n • uu 

22 


B 


= e'  E n,.'*®  • r‘ 
^ 1 la  la 


0) 


B 


(C-46) 


• g 

When  we  multiply  D_  , of  Equation  (C-44)  into  £ we.  get  two 

do  ~g  “ T>C 

terms  but  the  first  term  is  zero  because  is  an  n * n 

la  i-di 

skew- symmetric  matrix  of  vectors,  as  can  be  seen  from  making 
use  of  Equation  (B-17)  which  yields 

'I'B  . _ _jB  . 'I'B^  (C-47) 

ia  ia  ia  ia 


Evidently  the  matrix  on  the  left  is  equal  to  the  negative  of  its 
transpose.  From  this  it  follows  that 


(C-48) 


* B 

Hence  ^23^  a 


xa  ^ a 


m.?.®  • (C-49) 

^ a ^ 1 la  la 


Here  we  have  made  use  of  Equation  (B-19)  to  yield 


la  a a la 


(C-50) 


Equations  (C-46)  and  (C-49)  now  show  explicitly  that  Equation(4-81) 
is  satisfied ; i . e . , 


B 


0 


(C-51) 


-Cll- 


Evaluation  of  » 


From  Equation  (C-3)  and  (C-27)  we  immediately  get 

^ m.$.^  = /V  \ 

a 1 la  I a I ^a  a 

^ \^a  / 


where  we  have  made  use  of  Equation  (C-33)  in  the  form 


(C-52) 


= Z n,. 
) a 1 


la 
^ a 


(C-53) 


If  is  symmetric  (i.  e.  if  R.  and  its  first  two  partial 

la  la 

t 

derivatives  are  continuous)  then  7 t,  ® is  also  symmetric. 

pG  a 

^ a 

Evaluation  of 

From  Equation  (C-5)  and  (C-27)  we  get 

^ ^ c ^a  -I 

(C-54) 

O 

In  order  to  separate  u)  from  the  particle  dependent  factors  in 
the  second  term  on  the  right,  we  use  Equation  (B-15)  in  the  form 


uj®  • = uJ®  • 

ia  la 

Consequently,  written  as 


(C-55) 


=5^m.R.  • (f  • 2 rn  T®  • R^ 

32  “ 1 la  la  ■“  i la  1V.B  la  / a i la 

1 1 ^ ^ 


(C-56) 


-(  12- 


Evaluation  of 


For  0^2  obtain  from  Equations  (C-7)  and  (C-Z7) 


D33  • *fa 


' 4B 
.B  la  / ^ a 


$ 


la  1 


;B.  aB  . .B 


la  1 ia  ia 

1 


t,t  gi 


1 L ^ . 


$B  +IJB 

ia  ^ 1 la  la 


t 


(C-57) 


where  we  have  again  used  Equation  (C-55)  to  separate  from 
the  particle  dependent  factor. 




j Row  of  D^2 

Note  that  (v  is  a column  matrix  of  vectors,  and 

' ^ B la  / a 

^ a _ 

D is  also  a column  matrix  of  vectors,  whereas  D,,  is  an  n * n 
32  , ^ ^ 

matrix  of  scalars.  It  is  instructive  to  examine  the  j vector  ol 

D and  the  row  of  D,,.  In  order  to  do  this  we  note  from 
32  gt  t / r 

Equation  (C-22)  that  the  j row  of  (v  ^ ^ ) is  • 

^ T ^ a 


Thus  the  j vector  of 


lU 


B 


m . $ 


B 


i laj 


R,^ 

la 


(C-58) 


|.  

and  the  j row  of  is 

<D33>j. 


:B 

’ a 


la 


1 la 


. . 

laj  la 


a 

= e b'  V „ fc  i • ;b  . . , b' 

^a  i r.B  iajj  ia  ^ 1 laj  la 

1 ' ► 1 


(C-S9) 
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Note  that  by  writing  out  only  the  j row  of  (rather  than  the 

whole  matrix)  it  is  possible  to  separate  ^ from  the  particle 
dependent  factor. 


Evaluation  of  I 


According  to  Equation  (C-8)  is  equal  to  the  sum  of 
and  its  transpose.  ^22  given  Equation  (C-38)  and  therefore 

"En’iRia  ■ *fa  «a  ^ 


Now  note  that  the  third  terms  of  ^22  negatives 

of  each  other,  and  hence  these  terms  cancel.  Consequently 


T = +^m.R. 

a a 1 la  la  i i 

1 1 


~R  -R  ~R  ■**  ~R 

• C “ + ' 1 - I • 

la  ia  a a a 


Note  that  this  can  also  be  obtained  from  Equation  (A-58) 


in  the  form 


44  4^  ~R  ■**  ■**  ~B 

I =I  -I  '0U“ 

a a a a 


where 


n n 

\ - T.  il. 


- e'  V 

^ a a 

^ a 


and  where 
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V Jl  = Sni-  o R-  ^ R-‘  + Z-]  m.R.  • (7  r,  rL) 
a i\VB  lal  la  1 la  ' a-B  la/ 

t 1 E ^ 1 t _ 


= ’ R^  * m.R. 

1 la  la  I 1 


. 

la  la 


Evaluation  of  P ^ 
a 


(C-64) 


According  to  Equation  (C-9).  0 is  equal  to  the  sum  of 

— — t - ^ - 

^23  ■ ^3Z  given  in  Equation  (C-56).  is  given  in 

Equation  (C-44);  hence  is  given  by 

D,'=£n,.l‘>‘.  S'  =£m«B  • *®‘ 5=  + S'"  'i;  n,.  R,  • T® 

23  X la  la  ■“  i la  la  ^a  ” i la  la 


where  we  have  made  use  of 


= y^m.?.^  • r!^  + y^m.R.  • 

^ 1 la  la  ^ 1 la  u 


(C-65) 


Adding  i^23 

0^=2m.$®  •T.®£^+2m.R.  • ® • 0 ^ (C-bb) 

a ~ 1 la  la  a ~ i la  \ ,B  la  / a a 

1 ^ ^ 


plus  the  fact  that  • 0®  _ o;^  • (3  ^ . 

a a 

Note  that  Equation  (C-bb)  can  also  be  obtained  from 
Equation  (B-20)  in  the  form 

a a a 


CIS- 


whe  re 


L 


Bg  R 


^ j = i 

aj 


(C-69) 


B 


B 


The  transpose  ol  3 is  actually  more  easily  written,  and 


therefore  we  write 
n 


B„t 


B 


E 4 


B 


SfJ 


B 


. = b' 

a ’aj  j^B  ‘^.B^  a 


aj 


and 


Hence 


B 


B 


iyl) 


whe 


= En,.«B  • * E 


t t 


(C-70) 


= +Xm.R.  • V 

^B  a \*'^B  laj  la  ^ i ia  ^B  ia 


(C-71) 


(C-72) 


la 


(C-73) 


Bb 

Substituting  0 from  Equations  (C-72)  and  (C-73)  into 

“ . D 

Equation  (C-68)  yields  3^  as  in  Equation  (C-66). 
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J 


1 

I 


Evaluation  of  y 


, g 

According  to  Equation  (C-10),  y is  equal  to  the  sum  D,_ 
and  its  transpose.  ^33  given  in  Equation  (C-57),  and  therefore 

33  1 la  la  i la  ^a  la  i la  la 

1 1 L 4 - 1 

a 

(C-74) 

Now  note  that  the  second  terms  of  D33  and  negatives 

of  each  other  (because  of  Equation  (C-47)),  and  hence  these  terms 
cancel.  Consequently 

(c-75) 

a 1 \ la  / ^a  la  i la  ^a  ^.B  la  ' 

Is  1 s 

L^a  J »-  aJ 

We  can  also  get  this  directly  from  Equation  (C-13)  as  follows 


a V ® /a 


la 


• [?  5 


(C-76) 


We  hcavi-  now  obtained  explicit  expressions  for  all  the 
elements  of  D (see  Eciuation(  4-68  )or  (C-1)  ) and 'll  (see  Equation(  4-74) ), 
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